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Abstract 

The reformulation of field theory in which self-energy processes are 
no longer present ||] provides an adequate tool to transform 

Swinger-Dyson equations into a kinetic description outside any approx- 
imation scheme. Usual approaches in quantum electrodynamics (QED) 
are unable to cope with the mechanical momentum of the electron and 
replace it by the canonical momentum. The use of that unphysical 
momentum is responsible for the divergences that are removed by the 
renormalization procedure in the S'-matrix theory. The connection be- 
tween distribution functions in terms of the canonical and those in 
terms of the mechanical momentum is now provided by a dressing op- 
erator Ql that allows the elimination of the above divergences, as the 
first steps are illustrated here. 
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1 Introduction 



The obtention of the statistical mechanical properties of interacting quantum 
fields is the aim of studies in nonequilibrium quantum fields. Their relevance 
in describing the early stage of the universe in manifest but their interests 
lie beyond that application to provide an unified understanding of physics 
in all circonstancies. A kinetic description (of the Boltzmann-Langevin type 
of equations) requires the Swinger-Dyson equations as starting point from 
which a kinetic irreversible description is derived ||5|, via a generalization of 
the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy. Different 
schemes can be used to perform the transition from a reversible into an irre- 
versible description [^j. An appropriate truncation of the dynamics, through 
some molecular chaos assumption or a Kadanoff-Baym approximation, per- 
forms usually the job. In the framework of the scalar field, a adequate 
scheme of approximation enables to get rid of the divergences through re- 
summations (A^-points irreductible (N-PI) approximations [0]). The way of 
removing the infinities found in S-matrix theory, through the renormalisa- 
tion program initiated by S. Tomonaga, J. Schwinger and R.P. Feynman, is 
of no direct application since the description in the kinetic approach lacks 
the asymptotic limit used in computing the collision process. 

Our extension |^ of the dynamics of correlations provide a framework 
in which irreversible equations are derived in an exact way, that contains, 
as particular case, the reversible equations: the new dynamics is an en- 
largement from the original one that is entirely constructed from it, without 
introducing new parameters for instance. That approach enables to get rid 
of the self-energy processes at the level of the equations: their effect on 
the dynamics is taken exactly into account into the kinetic operator. We 
have dissociation of the approximation scheme and irreversibility. A further 
truncation of the hierarchy and approximation schemes would then be in- 
troduced inside an irreversible description and would not generate it. The 
situation is much more satisfactory on a conceptual point of view. 

Equations that no longer allow self-energy processes have to contain the 
solution to the divergences met in quantum field theory. The problem of the 
infinities is thought to be inherited from classical physics: the electromag- 
netic field dragged along with the charge of a moving electron turns out to 
provide an infinite mass to the particle. That last point has been challenged 
recently Q using the single subdynamics approach Q, From an 
application to classical electrodynamics [^, using the Balescu-Poulain for- 
mulation |1C] dealing with the mechanical momentum of the charged par- 
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tides, the usual mass divergence has been shown to be removable through 
resummations of classes of diagrams. A relation with the renormalisation 
group is manifest. 

Considering the unexpected sucess in the classical case, a new investiga- 
tion of the quantum case looks mandatory: does a formulation of QED exist 
so that no (infinite) substractions are required to provide a finite theory? 
We are convinced that a positive answer can be provided to that question. 
Two new ingredients are inserted with respect to previous attempts. The 
first one is the new mathematical tool constituted by the single subdynamics 
approach, the second one is the reintroduction of the mechanical momentum 
of the electrons (in opposition to its canonical one's) as one of the basic ob- 
servables of the description. A theory that deals with physical observables 
should be naturally finite! The difference between the canonical and me- 
chanical momentums in QED distinguishes the cases of the quantum scalar 
and vectorial fields. 

The relativistic interacting Dirac and electromagnetic fields necessarily 
involve an arbitrary large numbers of particles, namely electrons, positrons 
and photons. It is not possible to restrict ourselves to sectors involving only 
a finite number of particles. Therefore, it can be argued that such a system 
resort to quantum statistical mechanics and that the adequate methods are 
those developed in a non equilibrium statistical mechanical context Q. The 
physical systems, classical or quantal, are described by Wigner-like reduced 
distribution functions, depending on the position and the momentum of the 
particles (plus obviously their nature and spin) . The form of the infinitesimal 
evolution generator is dependent of the classical or quantal nature of the 
system. An interpretation in terms of usual probabilities is not possible 
in the quantum case: the Wigner distribution functions are not necessarily 
positive everywhere. 

We intend to construct a similar formalism in the high energy domain, 
using quantized fields. Reduced entities, involving a finite number of de- 
grees of freedom, are easily introduced by considering mean values of nor- 
mally ordered creation and destruction operators , or by taking the trace of 
the density operator on which destruction operators have acted on the left 
and creation operators have acted on the righttl. That procedure defines a 
complete set of moment functions that allows the computation of the mean 
value of any observable. The evolution equation of the set looks like the 

^Usual approaches use related n-points functions, introducing differences between the 
non relativistic and relativistic descriptions 
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B.B.G.K.Y. hierarchy. Since the creation and destruction operators of the 
Dirac field involve the canonical momentum, that set is defined in terms 
of that momentum. The Weyl's correspondence rule |^] between classical 
and quantal observables can be then used when the continuous observables 
are the position and the canonical momentum of the particles. Reference 
functions leading to a Wigner-like description are then introduced. Their 
argument involves naturally the canonical momentum of the particles. They 
constitute also a complete set for the description of the system. Their evo- 
lution equations can be computed from the evolution equation of the set 
of moment functions. The direct application of the single subdynamics ap- 
proach to the set of reference functions presents a priori no difficulty. When 
no external e.m. fied is present, a subset of reference functions, involving 
only the Dirac particle, satisfies closed equations of evolutions: the refer- 
ence functions functions involving the self-field do no longer appear in the 
dynamics. The analysis of the kinetic operator shows that it leads for the 
reference functions to the usual (logarithmically divergent) mass correction 
at the second order in the charge. 

However, our aim is a description in terms of the mechanical momen- 
tum. In Lorentz gauge, the form of the energy-momentum tensor of the 
electron-positron field shows that it involves also the electromagnetic field 
operators [11 1. That property can also be noted in the manifestly gauuge 
invariant formulation of quantum electrodynamics ||l^ where the distingo 
between the canonical and mechanical momentums is made explicit, unlike 
in more recent rediscoveries |13|. Moreover, the association between classi- 
cal functions of positions and mechanical momentums and the corresponding 
observables in second quantization is completely unknown and call for an 
extensive study by physicists. Two distinct roads are possible to obtain a 
description by Wigner functions defined in terms of the mechanical momen- 
tums. In the first one, a proposal for the above association is surmised from 
various considerations and the consequences of that choice examined. In 
the second one, it can be argued that the searched Wigner functions can be 
obtained from the reference functions since the latters form a complete basis 
of the description. The link is provided by a dressing operator acting on the 
reference functions Q. Inside the single subdynamics approach, when no 
external e.m. fied is present, a similar subset of the Wigner functions satisfy 
also closed equations of evolution governed by a dressed kinetic operator. 
It has been shown that, outside resonances processes a second order 
dressing operator modifies the kinetic operator at the same second order. 
Among the requirements for its determination, it can be required that the 
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new kinetic operator is finite. Therefore, we have to deal with a description 
the elements of which have a clear physical meaning, for true observables 
(e.g. the mechanical momentum), and so that these elements satisfy close 
equations of motion. The set of considered observables have thus to be com- 
plete and we have to scope with all physical situations. We can then expect 
that all physical entities are connected together by evolution equations free 
of any divergence. The resulting formalism is therefore rather heavy. 

In order to prove the feasability of the previous scheme, we consider in 
this paper a naive proposal for the above association between classical and 
quantal observables. We introduce somewhat arbitrarily a generalisation of 
the Weyl's correspondence rule. We are aware that high energy processes 
are not correctly included into that association. Therefore, our presenta- 
tion aims to accomodate possible alternative proposals and we focus on the 
structure of the new approach. It will be illustrated how the proposal de- 
termines a dressing operator that leads, for the second order corrections to 
the kinetic operator, to a logarithmic divergence analog to the second order 
mass correction, although no complete cancellation occurs. 

The usual description of the system in terms of a state vector in some 
Fock space is thus unappropriate for our purpose and, in the second section, 
we treat the electromagnetic field quantized in Lorentz gauge interacting 
with a (time dependent) classical source in the density operator formalism. 
That stage allows the introduction of all the required notations. Already at 
that level, the concept of superoperators (acting on the density operator) in 
normal form emerges. The Lorentz condition appears as a strong condition 
expressed in terms of the appropriate superoperator in normal form. Its 
persistence in time is explicitly shown. The properties of all classically 
generated incident electromagnetic field are thus made explicit. 

In §3, for the same system (classical source plus e.m. field), correla- 
tion functions are introduced for the description of the electromagnetic fied. 
They appear as the natural extension, in the density operator formalism, 
of the mean values of normally ordered creation and destruction operators. 
These correlation functions obey the equations of a B.B.G.K.Y.-like hierar- 
chy and the solution of those equations is provided: for a classical source, the 
correlation functions simply factorize. The non-classical aspects of the field 
have to be due to the non-classical character of the sources. Moreover, we 
show that the computation of the electric and magnetic fields mean values 
does not require the correlation functions involving the scalar or the longi- 
tudinal photon. A description similar to the quantization in Coulomb gauge 
is recovered: the fields depend on correlation functions involving transverse 
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photons and on the instantaneous Coulomb potential due to the classical 
source. 

The fully reduced description for the Dirac and electromagnetic fields is 
introduced in §4. The reduction bears not only on the electromagnetic fields 
variables but also on those of the Dirac field. 

In the next section §5, observables relative to one electron are consid- 
ered. The Weyl's rule of correspondence for quantum observables cor- 
responding to classical dynamical functions of the position and canonical 
momentum leads, from the consideration of their mean value, to our refer- 
ence functions expressed in terms of our previously introduced correlation 
functions. Wigner functions dealing for a complete set of particles observ- 
ables, including the mechanical momentum, are formally introduced. The 
energy-momentum tensor of the electron-positron field involves the electro- 
magnetic field operators and therefore, classical functions of the position and 
the mechanical momentum of the particle have no unambiguously correspon- 
dence with (in second quantization) quantum operators. That property does 
not prevent to present a simple (somewhat naive) association provided by 
a generalisation of the Weyl's correspondence rule. We are aware that high 
energy processes are not correctly included into that association that allows 
a first exploration of the properties of the corrrespondence. We think that 
the requirement of the finiteness of the final description can serve as a guide 
to provide that association. 

The one particle reference distribution function is further analysed in §6. 

In §7, the evolution equations for the whole set of correlation functions 
is derived. The number of elementary processes is of course much more 
important as in a non-relativistic case but all the contributions have been 
given an explicit form. 

The subdynamics approach is developped in the next section §8, includ- 
ing the introduction of the dressing operator. The equivalence of the single 
subdynamics and the usual approaches should not prevent unexpected re- 
sults using the new tool. The divergences met in the usual approaches may 
result from the use of an unappropriate description and we investigate here 
the first stages of that possibility. A finite kinetic evolution operator does 
not seem out of reach. Moreover, the link between non diagonal elements 
of the kinetic operators and the elements of the The dressing operatorpro- 
vides a conceptually satisfactory way of incorporating all the assets of the 
renormalization program of QED. 

Some conclusions are presented in the last section §9. 
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2 The density operator description for the elec- 
tromagnetic field 

We are concerned in this section with the electromagnetic field quantized in 
Lorentz gauge in interaction with classical sources, i.e. sources the evolution 
of which is a given function of time. This problem is well known in Hilbert 
space formalism but it is useful to reconsider it in the density operator 
formalism. 

Considering the undetermined number of field particles involved in the 
system, we use a statistical approach and the state of the electromagnetic 
field will be described in the next section by correlation functions, the evo- 
lution of which is provided by a BBGKY-like hierarchy of equations. 

The hamiltonian describing such a situation is well-known in the second 
quantization formalism we need [p^1- |ll9|| (see in particular the complement 
Ay in [^] from which the notations are borrowed): 

H = Hr + Hie (2.1) 

Hr is the field hamiltonian, decomposed in transverse {e,e'), longitudinal 
(/) and scalar (s) modes: 

Hr = Jd'k J2 ?i^kK(k)aA(k) + i)+y'd3fenu;k(a+(k)a,(k) + i) (2.2) 

X=e,e',l 

The interaction hamiltonian Hje is: 

Hie = ydV[-je(r,t).A(r)+cpe(r,t)^(r)] (2.3) 

The (external) charge density pe{r,t) and courant je(r, t) are here given 
functions of space and time that satisfy the charge conservation law: 

^/9e(r,t)+Vr.je(r,t) =0 (2.4) 

They transform as a quadrivector under a Lorentz transformation. 

The quantized potential operators A(r) and As{r) admit a plane wave 
expansion. 



A(r) = / + e-"<-4(k)) (2.5) 
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As is not hermitian with respect to the usual scalar product but with respect 
to the scalar product used in defining the (indefinite) metric. The infinite 
volume limit is used for the normalization of the creation and destruction 
operators of photons, characterized by the continuous wave vector k and the 
polarization A: 

[a,(k),a+(k')]=5^I,5(k-k') (2.7) 

K(k),a+(k')]=5(k-k') (2.8) 

The covariant quantization has been performed by using an indefinite 
metric, determined by a linear (unitary, idempotent) operator Af so that: 

MaxM = ax 

MasM = -as (2.9) 

The scalar product of two vectors \^ > and 10 > is then defined in the in- 
definite metric as < (plMlip > and the adjoint A of an operator A through: 
< (plMAltp >=< Tp\MA\(j) >*. Only the adjoint of the creation and destruc- 
tion operators for the scalar boson are modified by the indefinite metric and 
we have the link of the new adjoint (noted with an upper bar) with the usual 
adjoint (noted with a +): = —af and af = —as- Therefore, the scalar 
potential is, as previously noted, self-adjoint with respect to the indefinite 
metric. 

The mean value of an observable, hermitian with respect to the metric 
M, is defined, when the system is in a state {ip > hy < A >=< ^IM^I^/" > 
The class of admissible states is the class of states for which the Lorentz 
condition, that cannot be satisfied in a strong sense, holds in average: i 

< ^|J\M (^Vr.A(r) + ^i,^ IV' >= (2.10) 

The operator Ag corresponding to the time derivative of Ag, is defined in 
Schrodinger representation by: 

As = ^[As,H] (2.11) 
in 



^In reference |l^], that condition is erroneously expressed in the complement Av (58) 
by a mean value without the indefinite metric. The claimed equivalence between (67) 
on one hand and (66a) and the adjoint (in some sense) relation (66b ) on the other hand 



shows clearly that the Lorentz condition has to hold under the form (2.1C) 
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In a statistical approach, the system is described by a density operator p 
that has to be hermitian with respect to the metric. The factorized density 
operator pfac corresponding to a pure state > \s given by: 



and the mean value of an operator A is given by the trace in the usual 
expression: 

<A>= TV(A/j) (2.13) 

The hamiltonian operator is thus entirely determined in the second quan- 
tization formalism and, as pn, is hermitian with respect to the indefinite 
metric. The Liouville-von Neuman equation determines the evolution of 
the density matrix p/j for the radiation: 

ihdtPR = Lpn = [Hr + Hie, Pr]- (2.14) 

The density operator pR is constructed from factorized pfac density operator 
based on physical states satifying the condition ( p. 10 ). Superoperators allow 



a direct expression of that condition in the density operator formalism. They 
are defined as operators acting on the density operator. The Liouville-von 
Neumann (super)operator L ( p. 14 ) is an example. With an operator A 



acting on the left and an operator B acting on the right of the density 



operator, we define a factorizable superoperator A x B [20| by: 



{AxB)p = ApB (2.15) 



In second quantization, to each hermitian operator A, it can be associated 
pl| , Q in a unique way a superoperator A, written as a sum of factorizable 
superoperators in normal form (only destruction operators act on the left 
and only creation operators act on the right.) so that: 

<A>= Tr(A/9) = Tr(^p) (2.16) 

For instance, to the scalar potential Ag, we associate the normal form su- 
peroperator As defined as: 




(/ X a,(k))) 



(2.17) 
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and the superoperator As corresponding to Ag (2.11) is constructed in a 
similar way. The Lorentz condition can be expressed in a strong sense using 
the frequency positive components of the quadrivector A. In the present 
formaUsm, we have: 



1 



Vr.Air) + -As PR{t) = 



(2.18) 



The condition ( |2.10D i s recovered by taking the trace and using (2.12) for p. 
From ( 2.11 ) and (|2.3| ), we see that the density operator has to satisfy the 
following conditions for all wave numbers k and all times: 



(ai(k)-a,(k) + A(k,i))p^j(t) =0 
PR{t) (aKk)-a,(k) + A*(k,t)) =0 



where: 



A(k,t) 



-/9e(k,t) 



' 2eo?i(27r)3a;^' 

and Pe(k, t) is the Fourier transform of the charge density, i 



(2.19) 
(2.20) 

(2.21) 



Pe(r,t) 



(27r)3 



(i=^A:e*-''pe(k,t) 
Pe(k,t)= / d-'re-'^-^ pe{r,t) 



(2.22) 



The conditions ( p.l9 - p30 ) are mutually adjoint. Considering the self-adjoint 
character of the density operator, they are equivalent. For a time indepen- 
dent problem, condition ( 2.19| ) is equivalent to condition (66a) in Ay of |14]. 
If the conditions ( [2l9[]23c| ) are satisfied at some arbitrary time t, they re- 
main satisfied for all times. Using the charge conservation (2^), the proof 
is straightforward. 



3 The correlation functions for the electromag- 
netic field 

In a context of statistical mechanics of large systems, the relevant objects for 
its description are not the particles distribution function or the density 

^For the Fourier transform, we use the same convention as in in place of the con- 
vention of O], hence a difference in some numerical factors. 
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operator for the whole system but the classical or Wigner reduced distri- 
bution functions involving a limited number of particles. Considering the 
possibility of an arbitrary large number of photons involved, we treat the 
electromagnetic field along the same spirit and we have to define the ap- 
propriate reduced description. Such a description can be provided in terms 



a reduced density operator so that the mean value of any observable 
can be evaluated for all observables through a vacuum expectation value. 
This required that a normal form superoperator be associated with all ob- 
servables, as already considered in the previous section. It is possible to 
replace the reduced operator by its matrix elements that can be computed 
directly from the density operator itself and this point of view is adopted 
here. The state of the electromagnetic field is therefore not described by 
the density operator itself but by its set of correlation distribution functions 
pl| , obviously symmetrical in each set of variables: 

fnm 

= Tra,.i(ki)aj.2(k2) . . . a^^{k^)pn{t)a'^,J'k'^) . . . a+, (k2)a+ (k'l) 

(3.1) 

where the index u takes any of the values s, I, e, or e'. Those correlation 
functions appear as mean values of particular superoperators, those in nor- 
mal form, and depend on the diagonal, as well as the off diagonal in photon 
numbers elements of p. In the particular case where the density operator is 
factorized (the so-called pure case), those correlation functions are nothing 
but mean values, in Fock space, of a normally ordered product of creation 
and destruction operators. The normalization of the density operator re- 
quires in particular that /oo = 1 for all times. The correlation functions are 
not independent since they have to represent physical states and have to 
satisfy the relations due to ( |2.19| ) and ( 12.201) : 

/n+im(kiz/i; . . . ; kni^n', k„+i/; k']^z^^; • . . ; k^f^, t) 
-/n+im(kii^i; . . .■,knUn;Ki+is;'k[u[; . . . ;k^f^,t) 
-FA(k„+i,t)/„m(kii/i; . . . ;k„i/„;k'iz^i; . . . ;k^f^,t) = (3.2) 



fnm+i{k-ii'i', . . . ; k„i/„; k]^z^]^; . . . ; k^f^; k^_|_]^/, t) 
+/nm+i(kiz^i; • . . ; k„f„; k'^z^[; . . . ; k^z^^; k^^;^s, t) 
+X*{k'^+i,t)fnm{kiui; . . . ;k„z^„;k'ii/[; . . .;k'^u'^,t) = (3.3) 
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The electric and magnetic field operators B can be computed directly from 
the expression (|2.5| ) of the potential vector operator and from the form ( |2.1| ) 
for the Hamiltonian that is required for the evaluation of A(r) = i[A(r), H]. 
Their mean values involves only the moments /io(kA,t) and /oi(kA,i). The 
electric field operator is decomposed into the part Ej^, originating from the 
transverse component A_L(r) of the potential operator A (for A = £,£'), the 
part E;, originating from the longitudinal component A||(r) of the potential 
operator A (for A = /) and a part E^, originating from the scalar potential. 
Using the relations (|3.2[) and (|3.3D , usual manipulations lead to: 



< E;(r) >t + < E,(r) >t= -VrC/(r, t) (3.4) 

where U{r,t) is a Coulomb potential (without retardation), i.e. associated 
with the charge distribution taken at the same time t: 

The Lorentz gauge therefore provides expressions similar to those obtained 
in Coulomb gauge: the fields are derived from the transverse part of the 
vector potential < A^(r) >t and of the Coulomb potential U{r,t). 

Let us now consider the time evolution of the correlation functions. Prom 
the Liouville-von Neuman equation ( |2.14| ), we obtain immediately: 

^ /nm (ki 1^1 ;k2Z^2; • • • ; k„ zy„ ; k'^ i/J ; k'a ; • • • ;k^f^,t) 

= ^Trai,i(ki)a,,2(k2) . . . a,,„(kn)[HR + Hic,pR(t)] 
in 

xa+ (k:„)...a+(k'2)a+(k'i) (3.6) 

7n 2 1 

We have to substitute in this expression the hamiltonians and Hj^ by 
their expressions ( |2.2| ) and ( |2.3| ). It is useful to express the interaction 
hamiltonian Hje in terms of the Fourier variables: 




Hie = - j d^k \l 2eo{2TrfoJi^ {iei-K t).ei,x QA(k) + je(k, t).ek;^ a^(k) 



+c J d^k^ 2go(2t)3L^k ""'^^^ ~ ""^^^^^ ^^■'^^ 

In order to obtain a symmetrical expression, let us consider that the source 
and current densities define a vector with four components, labeled as the 
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photons: 

jsiKt) = js+ik,t) = jsAKt) = -CpeiKt) (3.8) 

We have to introduce a symbol s{iy) which takes into account the sign differ- 
ence for the scalar photon between the creation and destruction operators 
in (p.Tj): s(i^) is -1 if u = s and +1 in the three other cases. We have then: 



Hie 



2eo(27r)3cuk 

X (j,_(-k,t)a,(k) + s(z.)>+(k,t)a+(k)) (3.9) 

The use of the cyclic invariance of the trace leads then to: 
d 

^ /nm (ki 1^1 ;k2Z^2; • • • ; k„z^„; k^z^^ ; kgZ^a ; • • ■^'^m'^'m^t) 

i=l,n 4=1, m 

x/nm(klZ^i;k2Z^2; • • • \ knVnM\V \M2^'2\ ■ ■ ■ 'MrrAn^'^) 



V 2eo^(2vr)^a;k,^ 

x/n-im(kiz^i; . . . ; kj_iz^j_i; kj+it'j+i; . . . ; k„t'„; k'lZ^J; . . . ; kj„f^, t) 



-I 

i=l,m 



x/nm-l(klZ^i; . . . ; ^nt^n] ^l^^'l'i • • • j ^i+Wi+l'^ • • • j ^^f^, t) 

(3.10) 

This system of equations is equivalent to the Liouville-von Neumann equa- 
tion for the density operator. It can be solved in a simple way by recurrence. 
Indeed, in the r.h.s. of ( |3.10| ), only correlation functions with indices equal 
or lower to n and m appear. The equation for the first function /oo is trivial 
since the r.h.s. of ( |3.10| ) vanishes. That property corresponds to the trace 
conservation. We now consider the next functions /lo and /oi. Let us as- 
sume that the claasical charge and currents jy vanish before the initial time 
to- The solution of ( |3.10| ) can then be written in a trivial way: 



fu){^^lUl,t) =is{ui)J——^-- [ dne ^'^^i^* ^^^j (ki,ri) 
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/oi(k;z.;,t) = l \ - Tdrie^'K^* ^'^j,,{-k[,n) (3.11) 



Using the charge conservation law ( |2.4| ), a little algebra shows that the 
condition ( |3.2D is satisfied by our solution ( ^.11 ). 



The element /20 is now considered. The equation ( 3.1C| ) becomes for 
that element: 

d 

^/2o(kiz^i;k2Z^2,t) = -i (u^ki + ^^ka) /2o(kiz^i; k2Z^2, t) 
y 2eoft(27r)'^a;ki 



y 2eo?i(27r)'^a;k 



(3.12) 



The expression ( |3.11| ) for fiQ can now be inserted into that equation and it 
can easily be seen that the solution of ( p.l2D takes the factorized form: 

/2o(kiz^i;k2Z^2,t) = /io(kii^i,t)/io(k2J^2,i) (3.13) 

Such a factorization property holds for all the elements of / in presence of 
a classical source and we have: 

ui; k2Z^2; • • • ; k„z^„; k'li/i; k2i/2; • • • ; k^i^^, t) 
n hoiK^^,t)] i n /oi(k^^^.',t)) (3.14) 

The state of the field characterized by such factorization relations is called 
semi-classical or coherent. If an observable, associated with the emitted 
fields at various points, is constructed by the product of the (commuting) 
normal form superoperators associated with each component of the field, 
that factorization prevents a correlation betwen the emitted fields. The 
fluctuations in the emitted field are linked to the non-classical features of 
the sources. 

Let us examine more closely the normal form superoperators £(r) and 
B(r) associated with the field operators E(r) and B(r). For those super- 
operators, the same decomposition into scalar, longitudinal and transverse 
components is introduced, in analogy with the operators themselves . Ob- 
viously, the magnetic field superoperator B(r) only contains a transverse 
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contribution. The electric field superoperator £{t) contains in addition con- 
tributions arising from the scalar and longitudinal components. The sum of 
these two components can be given an explicit form: 



£.(r) +£,{r) =icj d'kJ ^J' v (e'k.r - a.(k)) x /] 



2eoa;k(27r)2 

/x(a+(k) + a+(k))l) (3.15) 



Considering the conditions (2.23-24), the density operator is an eigenopera- 
tor of £s(r) -|-£/(r) at all time: 



2eou;k(27r)3 

X (e*-'"A(k, t) + e-*^-'"A(-k, t)) pR{t) (3.16) 

Using the result ( |3.4D , the eigenvalue can be given a nice form: 

{£s{v) +ei{r)) pR{t) = -VrU{r,t)pR{t) (3.17) 



where the instantaneous Coulomb potential U{r,t) is given in ( |3.5D that is 
a particular case of ( [3.17 ). Therefore, as long as we restrict the observables 
to the fields, the scalar and longitudinal photons can be eliminated from the 
description and their effect taken into account by the instantaneous Coulomb 
potential. This proves the equivalence of the quantizations in Lorentz and 
Coulomb gauge. Such a property can also be extended to the case of multiple 
time measurements. 



4 The reduced description for the Dirac and elec- 
tromagnetic fields 

The approach of the previous section is extended to interacting Dirac and 
electromagnetic fields. The relativistic hamiltonian describing such a situa- 



tion is well-known in the second quantization formalism we need (|]ri|-|19], 
in particular the complement Ay (48-51) in [0] from which the notations 
are borrowed. 

H = Hd + Hr + Hj (4.1) 
where is the Dirac hamiltonian for the electron-positron field: 
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Cp and Cp are respectively the creation and destruction operators associated 
with the electrons characterized by the quantum number p that represents 
the canonical momentum vector p and the spin component a. The creation 
and destruction operators and bp are associated with the positrons, with 
a similar sense for the indices. Ep is the relativistic energy of the free particle 
(characterized by a bare mass m) in an eigenstate of the momentum (Ep = 
\Jm?'c^ A-p^c^). Hpi is the electromagnetic field hamiltonian, introduced in 
previous section ( |2.2p . The interaction hamiltonian Hj is the sum of the 
interaction hamiltonian Hje with external sources, considered in the two 
previous sections, and of the interaction hamiltonian -ff/D, involving (in 
normal form) the operators associated with the electron-positron field: 

Hid = J d\ [-jz)(r).A(r) + cpD(r)^(r)] (4.3) 

The Dirac field charge density operator pd{^) and current jD(r) are now 
operators given in second quantization (in terms of the Dirac spinor ^, its 
adjoint -0^ and the Dirac matrices oe.) {ip~^ is a spinor, the elements of which 
are the adjoint of the elements of tp, and its transposed is the adjoint of 

poir) = ge^[^+(r)V(r) - V^(r)V'+(r)] (4.4) 

jnir) = geC^['i/'+(r)aV(r) - iA(r)a V'*'(r)] (4.5) 

wher Qe is the electronic charge. The quantized Dirac field spinor ^(r) and 
its adjoint ip~^{r) are given by: 

V'(r) = [cpUpi^) + bpVpir) 
p 

4>+{r) = J2 [«(r) + bpv±{r)] (4.6) 

p 

The Up{r) are the positive energy plane waves spinors, eigenstates of the 
momentum and spin (along the momentum p) operators and of the Dirac 
Hamiltonian TCd for the free electron. The index p represents the set of 
quantum numbers (p, p) (p is the helicity index). The index p represents 
the set of quantum numbers (-p, -fi) that are opposed to those of p. The wave 
function Vp{r) represents a negative energy state and is therefore associated 
with the positron (to which is associated a positive energy). 
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The Dirac matrices are chosen to be given in terms of the Pauh spin 
matrices by: 



a - 



/3 







(4.7) 

The S and b matrices are given for future reference. We will also need the 
16 hermitian matrices that form a basis. If we define ctq as the unity matrix, 
that basis can be represented by B^ji = Oa® where ® represents the 
direct product of the Pauli matrices. In these notations, the component i of 
the vectorial or matrix is B\i. 

Fot the hamiltonians ( [4.1[]4.2| ), we use a continuum basis. The anticom- 
mutators of the creation and destruction operators are then expressed in 
terms of the Dirac delta functions. The quantized Dirac field takes then the 
form, taking into account the normalization of np^^(r) and t'p^^(r): 



V^(r) 



ft3(27r) 



d^v c^(p)np,^(r) + 6+(p)i;_p,_^(r) (4.8) 



Those operators c, c^, h and 6+ obey the usual fermion anticommutation 
relations and we have for instance: 



C/.(p),cJ(p') 



P - P 



(4.9) 



The explicit form of the normalized spinors we use can be given. The 
canonical basis for the spinors is constructed from the common eigenvectors 
to /9 and E^. Those vectors are noted |A/i > where A corresponds to the 
eigenvalue of /3 and /i to that of Ez. The link with the states in (|4.8| ) is 
|11 >= -uo+i, |1 - 1 >= uo-i, I - 11 >= t^o+i, 1-1-1 >= fo-i- The 
plane waves spinors Up^^j, Wp^^ can be obtained by the action of a unitary 
transformation T(p) |14] on those basis spinors: 

^^P,±i = r(p)no±i wp.ii = T(p)uo±i (4.10) 
where the unitary matrix is (p represents an unit vector aligned along p.): 

zp.r 



T(P) 



cos - (/^.p) sin 



exp ■ 



(4.11) 



where the angle Op is given by = Arctg(^). In terms of the relativistic 
energy Ep and the bare electronic mass, we have: 



6„ 



cos ■ 



2 En 



' E„ 



sm ■ 



2 En 



(4.12) 
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The r independent part of T is the unitary matrix t(p) (t^(p) = t(— p)): 



r(p) 



cos — (iSa.p) sin 



(4.13) 



The interaction hamiltonian can still be written in Fourier variables as 
with new expressions for the charge and current density: 

jx+ik,t) = (je(k,t) +jD(k,t)) .ek;, 
jx-iKt) = (je(k,t) +jD(k,t)).e_kA 

js {k,t)= 3s+ (k, t) = (k, t) = -c (/)e (k, t ) + PD (k, t) ) (4.14) 

The explicit form of the charge an current operators jDv± ( |4.14| ) in Fourier 
variables can now be given. 

jDu+{^) =11 j («J'+MM'(P'k)c+(p - ^;ik)c^/(p + ^hk) 



MM 



+a,"liuu'(p,k)c+(p 
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p, k)bf, 



hk)b^,{-p-^hk) 

1 1 
-P + 2^'^)V(P + 2^^) 



(4.15) 



and similar expressions for a minus subscript in place of a +. The detailed 
expressions of the coefficients a are provided in the appendix A. The charge 
and current operators are known |14| to satisfy the charge conservation 
equation that implies a link between the scalar density and the longitudinal 
part of the current. 

The hamiltonian operator is thus entirely determined in the second quan- 
tization formalism and we can now consider the equation for the density 
matrix p that satisfies the Liouville-von Neuman equation: 



ihdtp = Lp= [Hd + Hr + Hi,p] 



(4.16) 



The density operator p has to involve only physical states. Therefore, it has 
to satisfy the following conditions (2.19- 2!2C| ) for all wave number k (c.f. [^] 
(67) p. 425) where A(k, t) is now a self-adjoint operator in the Dirac field: 



A(k,t) 



(pe(k,t)+pD(k)) 



(4.17) 
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where PdO^) is now the Fourier transform of the charge density operator. 
As for external sources, if these conditions ( 2.1£ - 2l20| ) hold at some given 
time, they hold for all times as a consequence of the charge conservation. 

The basic entities for our description of the interacting electron-positron 
and electromagnetic fields are defined, as in (3.1), by correlation functions 
obtained by mean values of particular superoperators in normal form: 



Jll...^ls^l'|...^l[jll...jl-sjl',...^^\. , , , _ _ 

f ss'ss'nn' (pi ; P2 • • • Ps ; Pi ; P2 • • • Ps ; Pi ; P2 • • • Ps ; 

pi; P2 • • • P^; kiz/i; k2i^2; • • • ; k„z^„; kii/[; k'ai/a; • • • ; ^n'^'n'\t) 
= ;i-i(^+^'+^"+^-')Tra,,(ki)a,,(k2) . . . a,„(k„) 



x6^^(pi) . . .6^-(pg)c^^(pi) 



{Vs)p{t)ct {^',, 



■<i(Pi) 



6+(p;)a+ (k:,,)---<(k'2)<;(k;) 



(4.18) 



The variables with an upper index "tilde" are associated with the positrons. 
The variables with a "prime" are associated with the creation operators, as 
in (3T). The element /oooooo corresponds to the normalization of the density 
operator and takes the value 1. When the arguments of the f ss'ss'nn' func- 
tions are labelled as in the definition ( [4.18| ), they will not written explicitly. 

It is useful to introduce substitution operators S for the variables of 
electromagnetic or electron-positron fields. Acting on a function of kj, z^j, 
the operator S'(k, i/; kj, fj) replaces the variables kj, fj by new variables 
k, i^, and a similar definition for the matter variables. Those substitution 
operators allow to avoid to write explicitly the variables of the f ss'ss'nn' 
functions: only the modifications with respect to those of the definition 
(|1|) have to be given in terms of the substitution operators. 

These moments ( [4.18| ) are not independent since they have to represent 
physical states and have to satisfy relations due to ( |2.19| ) and ( p. 20 ) besides 
those induced by quantum symmetry of the Dirac field From( 2!T9| ), the 
following relations have to hold: 



?i-5(^+^'+^"+^"')TVa,,(ki)a,,(k2) . . . a,„(kn)b^,(pi) . . .b^,(ps 



(pi)...c^ 



aiCk) - as(k) + 



(pe(k,t)+ :pD(k) :) 



Pit) 



XC+ (p^O . . . c+ {p[)b;,_ ip',,) ... 6+ {p[)a+ (K,) . . . a+ (k^)a+ {k[) 







(4.19) 



from which, using the cyclic invariance of the trace, relations for /, 
that generalize (|3.2[j3.3|) can be deduced. 



ss'ss'nn' 
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For a classical external source, the longitudinal and scalar components do 
not appear in the computation of the mean value of the electric field (|3.4| ). A 
similar property holds for the combination of quantum and classical sources. 
If we neglect for the moment the time dependence in the classical source, 
to the electric field observable, we can associate a superoperator in normal 
form. Its transverse components are given by analogy with the external 
source case. The sum of the longitudinal and transverse components is still 
given by ( 3.15| ). The use of the relations ( p. 19 -2. 20) in ( 3.15| ) imposes a 
reorganization to obtain an operator in normal form: 



1 



, k / 

(2vr)2 2eo7 u;^ V 



J (fp( all ^>{p,k) 



1 



1 



V(P + o^^) X c+(p - -hk) 



+asMM'(P''^) 



b^,{p + -hk)xb+{p--hk) 
Ixc;{p-Uk)b+{-p-^hk) 



1 



1 



p + -hk)cn'{p + -hk) X / 



(4.20) 



For all types of sources, the mean values of the electric field operator can 
therefore be easily expressed in terms of the correlation functions that in- 
volve neither the longitudinal nor the scalar photons. 



5 The observables of the Dirac field 

The continuous observables associated with the Dirac field are the positions, 
the canonical momentum and the mechanical momentum, that can be com- 
bined with spins and nature (electron-positron) properties. The mechanical 
momentum is more physically relevant than the canonical momentum but, 
for commodity reasons, the correlation functions ( 4.1^ ) are more easily de- 



fined in terms of the canonical momentum. For a physical interpretation, 
we have however to keep in mind that only the mechanical momentum is 
meaningful and that the canonical momentum has to be considered as a 
useful intermediate variable. 

In non relativistic quantum theory, to a classical function Odiic, p) of the 
position and momentum of the particle, is associated a quantum observable, 
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O defined through the operator (Weyl's correspondence rule) 

0=J^J rf=^Kd3xaKK,X)e'(K.f+iX.p) (5 -^) 

where Oci(K.,^) is the Fourier transform of Oc/(r, p): 

Oe/(K,X) = 1| dVd3pO,,(r,p)e-'(K.r+iX.p) ^g^) 

The second quantization form 0[2Q] of an operator O is constructed a priori 



as: 



[2Q] 



dVV+(r)OV(r) 



(5.3) 



The present relativistic approach requires the generahzation of those con- 
cepts. No straightforward procedure is directly available. Indeed, the ex- 
pression of the mechanical momentum, derived for instance from the energy 



momentum tensor [11 1 p419 requires also the fields operators describing cre- 
ation and destruction of particles, besides the operators dealing with the 
field. 

However, for the canonical momentum, no such problem arises. There- 
fore, we consider the Weyl's procedure for that momentum and possible 
extension to the mechanical momentum will be presented. 

The joint characteristic operator 'E^^\r,v,a, f3), involving the position 
and the canonical momentum of the electron, and the hermitian operator 
Bafs previously introduced (Bafj = (Tp), is defined through an expres- 
sion that respects the symmetry between electron and positron. Due to 
the presence of a general Ba/3 matrix, the simple substraction procedure 
present in (|4.4| ) , (|4.5| ) cannot be generalized and the normal ordering of the 
creation and destruction operators has to be restored by hand, without con- 
sidering the contribution due to the anticommutator. The existence of that 
procedure is noted by a prime on the bracket. For one particle, we have 
(The generalization to an arbitrary numbers of particle is staightforward 
and allows to introduce a complete set of observables.): 



H(o)(r,#,a,/?) = e*(K.f+iX.p)^^ 



a/3 



[2Q] 



(fr 



X.V+iK.r 



-Ba/3V'(r) 



(5.4) 



The Weyl's procedure cannot be applied directly in the case of the me- 
chanical momentum. A naive proposal based on it provides nevertheless 
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clues on the structure and potentialities of the approach. In a non relativis- 
tic approach, the mechanical momentum of the electron w is linked with the 
canonical's one p by the relation : 

v=p-qeAj_{r) (5.5) 

where the transverse vectorial operator A_|_(r) can be deduced from ( [2.51) . 
The use of ( |5.5| ) in a relativistic context in place of the true momentum 
tensor [0] p419 provides a limitation in specific results. 

A_L(r) contains creation and destruction operator associated with the 
field. The use of a similar expression in the Dirac description of the electron- 
positron particles requires some adaptation. Indeed, the description involves 
simultaneously both kinds of particles and we have to provide a mechanism 
that enables the correct sign of the associated charge that plans a role in 
( |5.5D . A free electron is represented by the positive energy plane waves spinor 
Up{r) {\up > in Dirac notation). We can therefore introduce a projector 
{Pe = J2p\up >< Up\) on these states and a similar projector Pp for the 
positrons. 

The basic observables for the Dirac field are the position and the mechan- 
ical momentum of the electron-positrons, and the hermitian operator 
introduced previously We have to propose a form for the joint characteristic 
operator involving these basis observables. 

The hermitian operator B^f^, is defined through an expression that re- 
spects the symmetry between electron and positron. The difference between 
the canonical and mechanical momentum involves the transverse field opera- 
tor A_|_(r) and the operators Pe and Pp. The ip{r) operator is not symmetri- 
cal with respect to the electron-positron symmetry. Indeed, the destruction 
operator associated with the electron would be on the right of the gener- 
alisation of (|5.3D , as the creation operator for the positron. In presence of 
a general i?^^ matrix, the normal ordering of the creation and destruction 
operators of the particle has to be restored by hand, without considering the 
contribution due to the contractions between similar fields. The existence 
of that procedure is noted by a prime on the bracket. We want a similar 
association for both kinds of particle. We separate the positive frequency 
part A|^^(r) (with the destruction operator) and negative frequency part 
a||_ ■*(r) of A_L(r). If A'|^^(r) has to be associated with the electron de- 
struction operator, a||_ ''(r) has to be associated with the positron creation 
operator. For the term linear in A^(r), it will be manifest that the field op- 
erator has to be outside the first quantisation operator e*''^'''"^^"^^'* where 
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p will stands for jV. The operator to be placed at the right will be noted 
A±r{r) and the operator to be placed at the left will be noted A±i(r). We 
have 

Ax.(r) = Ai+)(r)Pe - A'^^\r)Pp (5.6) 
An{r) = A^^\r)Pe - A^]^\r)Pp (5.7) 
The proposal for the joint characteristic operator H(f,#, a, /?) is: 

^,(K..+lX.*)^^^|^^^^ ^1^3, [^+(r)e-^^-^-«^X.V+.K.r^^^ 
xe-^X.A,„(r)^(r)]' (5_8) 

The choice of the correspondence has to be justified later on by the properties 
of the description that it generates. We anticipate by noting simply that 
if we interchange A_|_/(r) and A±r{r) in ( |5.8| ), no contribution analogous to 
the second order mass correction is obtained. 

As we have seen in quantum optics Q, a dressing operator can be in- 
troduced in the single subdynamics approach. That dressing operator will 
be able to connect all the descriptions arising from different proposals for 
the connection between classical and quantum observables. In the present 
study, we will simply consider that proposal ( |5.8| ) and examine mainly its 
consequences at the level of the finiteness of the evolution operator. Rel- 
ativistic corrections may be required to the present proposal based on the 
non relativistic Weyl's rules. For a proposal that is not based on the energy 
momentum tensor, two different gauges may provide inequivalent represen- 
tation: we have no reason to suspect a gauge invariance of our first proposal 
that is formulated without caring about that point but aims to ilustrate the 
feasibility of the approach . 

The transverse component (also the longitudinal one) of the field op- 
erator A_L(r) does not commute with powers of the operator V and that 
expression (^^) is defined through the perturbation expansion of the expo- 
nential. 

As we have seen in quantum optics 0, a dressing operator can be intro- 
duced in the single subdynamics approach. The dressing operator Q is a 
tool to perform the transition between the distribution functions in terms of 
positions and canonical momentums towards distribution functions in terms 
of positions and mechanical momentums. It is entirely determined by the 
proposal (|5.8|). That dressing operator is able to connect all the descriptions 
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arising from different proposals for the correspondence between classical and 
quantum observables. 

For the mechanical momentum #, the simple form ( |5.8| ) is not expected 
to hold in general and, for all proposals, we would obtain, for the joint 
characteristic operator H(f a, the following structure {ki is here a 
shorthand writing for the set of variables ki, z^i): 

H(r,#, a, /?) = ^ ^ y" d^p J J d^ki ... j d^kn J d^k[ ...J d^k'^, 

n,n' fifi' 

X E E (7^^''^'""'(p,p',A:i,...A:„;fcl...C)<,(p')c,(p) 

-L n' 

xa+ (k;,) . . . a+, (k^)a+, (k;)a^^(ki)ai.2(k2) . . . a^„{\<in) 

n' 2 1 

^^WM'nn'^p^ P', fcl, . . . k^- k[... k'^,)b+{p')b,ip) 

xoi (k;,) . . . a+, (k'2)a+, (k[)a^^{ki)a^,^{k2) . . . a^,„{^n) 

n' 2 1 

+7^^''^'""'(P, VM,... k^; k[... C)<(p)6j (pO 
xa+, (k^,) . . . a+, (k^)a+, (k;)a^^(ki)a^2(k2) . . . a^^^iK) 

n' 2 1 

xa+ (k;,) . . . a+ (k^)a+ (k'i)a,,(ki)a,,(k2) . . . a,„(k„) + ..\ (5.9) 



The dots refer to terms involving more particles creation and destruction 
operators. These 7 functions define the joint characteristic operator for 
one particle. Their dependence in K and X is implicit. The knowledge of 
these functions is still a challenging problem for physicists. Even in ordi- 
nary quantum mechanics, the correspondence between a function of classical 
observables, so that their quantum operators do not commute, and the asso- 
ciated function of quantum operators is not settled: Weyll correspondence 
rules is only one useful proposal. In the relativistic case, the mechanical 
momentum of one electron can be obtained throught the energy-momentum 
tensor but at our knowledge, no proposal for products of that momentum or 
the joint distribution of product of the moment and the position operator 
has ever be emitted. Imposing the finiteness of QED involves constraints to 
the 7 functions. We think that these constraints are natural and we link in 
fact the finiteness of the description with the general correspondence rules 
between a classical function and the associated quantum operator. Gauge 
invariance provides of course also constraints on the choice of the 7's. 
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To the form (5.4) correspond a limited number of non-vanishing 7 func- 



tions (both n and n' vanish and no supplementary terms are present) that 
can be easily computed. They will constitute our reference description. The 
form (|5.8| ) defines non trivial 7 functions and enables to illustrate the use 
and properties of (5.9). 



The mean value of the characteristic operator for the canonical momen- 
tum 'B^'^\r,w,a, (3) involves the trace of the product of that operator with 
the density operator and can be written in terms of the moments that we 
have previously introduced ( [4.18| ), with the same convention for the name 
for their argument: 



H(°)(r,#,«,/3) =;.3 / d'p J dV7r'^'°°(p,p')/r/o'ooo 
\fj.fj.' 

jlp.' 

+ E/ d'p' I d'p'jlf''-' {p',pV{;[Lo 

+ T.Jd'pJ d'p^lf-^'\p, mLoo\ (5-10) 

If we consider the quantum operator OBap associated with a classical func- 
tion Ocz(K,X) and the quantum operator B, its mean value is provided 
by 

< OBap >= J d'Kd''xdci{K,X) < S(i-,#,a,/3) > (5.11) 



where it represents the mechanical momentum of the electron-positron field. 

fW 

ial3 



We can define the generalization of the Wigner distribution functions 



by imposing 

< OB^p >= J d^Kd'xd,i{K,X)f^p{-K,-X) (5.12) 



and by identifying the two previous expressions ( 5.11 ) and ( ^.121) . The use 



of H^^)(r,ir, a, j3) in ( 5.11] ) and a similar identification provide the definition 



of the reference functions f^^- 
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The 16 dressed distribution functions are the natural variables to 
deal with the electron-positron field for observables involving only one par- 
ticle. They are completely determined by the gamma's and expressed in 
terms of the bare distribution functions fs^gs'nn' values of n and n'. 

If the observables involve more than one particle, we need also joint Wigner 
distribution functions that can be defined in a similar way. 

These relations ( |5.1C1| )-( |5.12| ) and the similar ones for more that one par- 
ticle (for describing for instance a collision process) form a basis for studying 
the time evolution of the joint dressed distribution functions describing 
all physical observables associated with the particles. 

Although we deal with a complete set of quantum observables for the 
particles, the set of Wigner distribution functions for the particles is not 
closed a priori by the evolution equation: indeed, distribution functions in- 
volving the em field appear naturally. When no incident field is involved, 
our reformulation of quantum field without self energy parts allows to ex- 
press those functions in terms of the distribution functions involving the 
particles only through a so-called creation operator. Therefore, in that ap- 
proach, closed kinetic equations can be obtained for the Wigner distribution 
functions when no incident field is needed. That incident field represents of 
course an external field and not the self-field generated by the particles: the 
effect of that self- field is already taken into account in the kinetic equations. 

Once the gamma's (5^) are known, the link provided by the creation 
operators can be used to express < H(f,#, a,/3) > in a way similar to the 
right hand side of ( 5.12| ): only distribution functions of the particles are 
present. They can be expressed in terms of the reference distribution func- 
tions /^^ describing the canonical momentum. That process is equivalent to 
the determination of the set of Wigner distribution functions in terms 

of the reference distribution functions /^^'' , relative to the canonical momen- 
tum, obtained in using ( ^.10 ) in place of ( ^.11| ). It corresponds thus to the 
determination of a dressing operator. We know Q that a dressing operator 
of the second order (in the charge) will modify the kinetic equations at the 
same order, outside resonances. We will see in a next section that the second 
order kinetic operator contains the usual divergent mass correction for the 
reference distribution functions fj^^. The supplementary contributions due 
to the dressing could compensate that divergence. Our attempt (|5.8|) that 
aims to replace exponential of the canonical momentum p by some combi- 
nation of exponentials of p and geA(r) in ( ^.Ij) does not provide the solution 
since high energy processes are not correctly described, due to the lack of 
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fermions pairs creation and destruction. That attemp proves nevertheless 
that new terms, presenting the same kind of divergences, appear although 
no complete cancellation occurs. 



6 The one particle reference distribution func- 
tion 



The explicit evaluation of the first terms in the joint characteristic operator 
(5.4) is the object of this section. In place of (5.4), we can use a more 
symmetric form with respect to the electrons-positrons field using (Vop is 
a derivative operator acting on the left cf. [0]). The final result do not 
depend on the form chosen 



H(?,ir,a,/?)|(°) = e^(K.f+iX.p)^^ 



a/3 



^+(r)eX-(|(V-V„,)+iK.r^^^^(^^ 



(6.1) 



We see that expression as the zero*'' order in the charge of H ( |5.9| ), hence 
the upper index. 

The well known relations: 



X.V+iK.r _ ilK.X iK.r^X.V 



(6.2) 



will be useful for the computation of ( |6.1[ ). Since [^(V— Vop), r] = [V, r] = I, 

and [V, Vop] = we have also 



X.i(V-Vop)+iK.r ^ giiK.X^*K.r^X.i(V-Vop) ^ ^-X.i(Vop)^iK.r^X.i(V) 



6 2 



'6 2 



(6.4) 



Using an integration by parts, we obtain at zero order in the charge: 



gi(K.r+iX.ir) ^ 



a/3 



(0) 
[2Q] 



V;+(r)6^K.r^^^^(^^X) (6.5) 



Using the expression ( ^ ) for ^/;(r) and the relations ( [4.11D , ( |4.1C| ) and ( [4.11D , 
we obtain 



V'(r) 



1 



ft3(2^) 



E / ^''P L(p)T(p)no^(r) + 6+(p)r(-p)^;o-;.(r)l (6.6) 
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We note r(p) given in ( 4.11| ) as 



r(p) = r(p) exp 



zp.r 



where the unitary matrix t(p), given by ( 4.13| ), takes the form 



cos ((g)/) — i sin Y (cj2 cr.p) 



(6.7) 



(6.8) 



and T^(p) = r(— p). We use also the previously introduced notation // 
—11. Therefore, 



V'(r) 
'0+(r) 



1 



1 



ft3(27r 



ip.r -ip.r 
Cf,{p)e fi r(p)uo;, + 6+(p)e « T-(-p)uo;i 



-ip.r ip.r 
c;r(p)e^^uoMr(-p) + bf,{p)e~vopT{p) 



(Pi P j IS given by 



(6.9) 
(6.10) 



(P, P ) 



jiK.X / ,3 



re 2 



-ip'.r 



tp.(r+X) 
xe -r(p)tio;. 



(6.11) 



The integration over r can be performed and provides a Dirac delta function. 
We introduce a notation for the geometrical factors: 



C^^(P ,P) =< eV|r+(p')i?a/3T(p)|e/. > 



(6.12) 



The links between the states |.^;U > on one hand and uq^, vq^ on the other 
hand has been provided previously. Therefore, 

ip.X 



llf'''\p,p') = ^e4K.X^(p-p' + ;.K)e^e;^"(P >P) (6-13) 
Similarly, 

lap'^ (p,p) = ^3- 



.Woo^p^p,) = -le^^K-^5(p-p' + nK)e^^f' (-P,-P') (6.14) 
.>^-'°°(p,pO = i,e^^K.X,(p + p, _ nK)e^5lf'{p, -p') (6.15) 
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lAtVOO, 



(-p',p) (6.16) 



The gamma's corresponding to the reference one particle distribution func- 
tion is thus explicitly known. They allow to express the reference functions 
in terms of the correlation functions ( 4.181) that do not involve the field. The 
two sets of functions provide an equivalent description. The gamma's corre- 
sponding to the proposal ( p.8| ) can be found in appendix B. The connection 
between the a's (4.15) and the (5's ( |6.12 ) can be made explicit. For v ^ s, 
we have the relevant values: 



a 



11 



(p,k) = qeC^ {ey)i5\f^'{p - \fik,Y> + ^ftk) 



a 



11 



(P,k) 



a 



11 



a 



11 



i=l,3 



i=l,3 



;p, k) = qec^ (P - ^^k, p + iftk) 

1=1,3 

:p,k) = gecp 5](e,),<5;fi^'(p-iftk,p + i?.k) (6.17) 



j=l,3 



The next section is devoted to the time evolution of the correlation func- 
tions. 



7 The evolution equations for the correlation func- 
tions 



The evolution equations for the correlation functions (4.18) can be obtained 
from the Liouville-von Neumann equation ( 4.16| ), using the cyclic invariance 
of the trace. The structure of the equation is the following: 



dt 



ss'ss'nn' 



< ss ss nn \L\ss_ssnn > /. 



' t? iV) '>-l / 



ss' s s nn 



(7.1) 



s s's s n n' 



The evolution operator C admits the same decomposition as L ( [4. 161) or H 
( [4.11 ). From the explicit form of the Dirac hamiltonian ( |4.2| ), written in the 
continuous basis as: 

^D = ^J2l d^E, (c+(p)c^(p) + 6+(p)6^(p)) (7.2) 
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it can easily seen that the contribution of Co is diagonal and its only non- 
vanishing elements are: 



< ss'ss'nn'\CD\ss'ss'nn' >= ^ \ ^2 ~ ^p' ~ 



ii=l,s i=l,s i=l,s' 



i=l,s' 







The contribution of Cr is also purely diagonal and corresponds to the first 
term of (3.10): 



< ss'ss'nn'\£D\ss'ss'nn' >= —i ^ a;^. — ^ a;^' ) C^-^) 

, j=l,n 



i=l,n' 



The contributions due to Cj are not diagonal. The non-vanishing contribu- 



tions arising from the external sources can be inferred from (3.10): 

-S{l'i)3eu,+ {^i,t) 



< ss'ss'nn'\Cie\ss'ss'n— In' >= i ^ 



1 



X I ]J S{]<ij+i,Uj+i\]<ij,Uj) 

\j=i,n-l 



(7.5) 



< ss ss nn \Lie\ss ss nn — 1 >= 




2eoft(2^)3c^k. 

xs{ui)jeu,-{-Kt) { Yl S{k'j^^,uj^^;k'j,uj)\ (7.6) 
Vi=«,"'-i / 

The contributions due to Cjd are more complicated since they involve 
changes both in the e.m. field and the Dirac field variables. In order to 
visualize those contributions, it is useful to introduce a diagrammatic repre- 
sentation for them. Let us represent the correlation function fss'ss'nn' by s 
left orientated straight lines, bearing a plus sign, s' right orientated straight 
lines, bearing a plus sign, s left orientated straight lines, bearing a minus 
sign, s' right orientated straight lines, bearing a minus sign, n left orien- 
tated wavy lines and n' right orientated wavy lines. The different elements 
of CiD can then be represented by vertices inducing changes in the corre- 
sponding lines. These vertices involve at most one photon. Each occupation 
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numbers s, s' , s and s' can change by at most one unit and the modifica- 
tions are restricted to a couple of them, leaving unchanged the parity of 
s + s' + s + s' . Moreover, the modifications leave invariant S, given by the 
expression S = s — s' — s + s'. S corresponds to the difference of the charge 
number A'^^: = s — s oi the state at the left of the density matrix with the 
charge number Nji = s — s of the state at the right. Each charge number 
is no longer individually conserved, as in the Hilbert space or in the density 
operator description, but only their difference is invariant in the correlation 
functions formalism. Each value of S defines therefore a sector and the evo- 
lution is completely decoupled between those sectors. The equation ( |7.1| ) 
can therefore be analysed in each subdynamics defined by those sectors. 

Those equations ( [7.1D , with the explicit values of the vertices (given in 
appendix C), together with the expression of the mean value of the physical 
observables ( [5.12|) , form the basis for a further analysis of the evolution of 
the coupled Dirac and electromagnetic fields. 



8 The single subdynamics approach 

8.1 The extension of dynamics and the choice of the subdy- 
namics 

We intend to describe a Dirac field in interaction with the electromagnetic 
field. The structure of our description is the same as in the classical case 
and we proceed in a very similar way. From our hierarchy equations 
of motions for the correlations functions ( [4.18 ), we derive a description in 



which we can distinguish whether a field is incident, emitted or dressing while 
including the original description. Of course, such an extension is largely 
arbitrary: the criterion may indeed be defined in various ways. Since we 
do not depart from equivalence conditions with the original formulation, 
all choices are equivalent and are simply determined by commodity reasons 
for the physical situation in mind. In the classical case, the distinction 
between incident, emitted or dressing field could easily be based in reference 
with each classical charge [^J . Considering the quantal undistinguishability, 
leading to relations induced by the quantum symmetry [^, it is convenient 
to define the criterion with respect to the whole set of Dirac field. Each 
photon is therefore labeled by a supplementary index rj that describes its 
status, either as incident photon rj = i, outgoing photon r/ = o or dressing 
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photon (self-interaction) rj = M. The self interaction includes here transfert 
of photons from a charge to "another" one. Since our main interest is to 
investigate the presence of infinities due to the self interaction, the above 
criterion is suitable for our purpose. 

The equation of motions for the new set of correlation functions 

f ss'ss'nn' (pi ; P2 • • • Ps ; Pi ; P2 • • • Ps ; Pi ; P2 • • • Ps; 

P'i;P2 • • •P's;klZ^l??i;k2Z^2??2; • • ■ ;Kl'nr]n;K^Wl'^'^W2V2^ ■ ■ ■\^n'^'n'^'n''^t) 

(8.1) 

is constructed in the usual way: the non-vanishing matrix elements are 

numerically the same as in fg^/g^nn'' ^^^^ certain matrix el- 

ements vanish and reflect the nature of the photon absorbed or emitted: 
for instance, an incident photon is never emitted by the interaction and a 
outgoing photon cannot be absorbed. The constitutive equations are similar 
to the classical case: 

Jss'ss'nn' (Pi;P2 • • •Ps;Pi;P2 • • •Ps;pi;P2 ■■■Ps; 

p'l; P2 • • • Ps; kiz^i; k2J^2; ■ • • ; Ki^n, kX; ^1^2; • • • ; ^n'^'n'it) 

Jss'ss'nn' (pi ; P2 • • • Ps ; Pi ; P2 • • • Ps ; Pi ; P2 • • • Ps ; 

{v} 

P'l; P2 • • • P's! kiz^ir/i; k2Z^2?/2; • • • ; k„i/„7?„; \l^v'^ri^\\l^v'^rl^\ . . . ; \l^,v'^,rl^,\l) 

(8.2) 

where the prime on the sum means that only the values of 77 =i or o have 

to be included. Providing that the f^lig^'Xi'' ^^^^ functions keep 

their value upon substitution of the value o of any rj by the value s (and 
vice-versa), the set of equations for / provide for /, via the constitutive 
relations, the original set of equations. Of course, the new set is richer since 
the link between the initial conditions is not yet determined. In the classical 
case [^, that link has been provided by the condition of belonging to the 
single sub dynamics and that condition holds also here. 

We have to define now our subdynamics and define the "vacuum" of 
correlations. As in previous papers, the vacuum is defined by all the "states" 



*The procedure is reminescent of the in and out states usual in S-matrix theory but no 
aymptotic limit is required here for the validity of the approach, hence the consideration 
of unstable states does not prevent any peculiar difficulty 
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that do not include a self-field: a / function belongs to the vacuum if all the 
indices rj are of the i or o type. The set of equations is not close a priori 
for the corresponding elements, if the criterion of belonging to the vacuum 
bears only to the self-field. Indeed, in the kinetic equation, the one particle 
Wigner distribution function is driven by the value of all possible Wigner 
distribution functions involving an arbitrary large number of Dirac particles. 
Indeed, all those functions are considered as independent and we deal with 
the equivalent of a BBGKY hierachy. 

Other choices may be more appropriate acccording to the physical sit- 
uation in mind, e.g. if the aim is not the description of an assembly of 
electrons-positrons, with an arbitrary initial condition, but the description 
of only a few specified Dirac particles, for which the initial condition for the 
Wigner distribution functions is known. Through another choice of the vac- 
uum, a solution of the hierarchy equation can be obtained so that it is driven 
by the elements describing those specified Dirac particles, while the other el- 
ements, also present, are determined as time-independent functionals of the 
relevant one's. They describe the "cloud" of virtual particles, determined 
by the "real" ones. For the elements of the subdynamics approach that will 
be considered here, the specific choice is in fact irrelevant. 

Let us present our scheme for the determination of the closed equations 
for the Wigner distribution functions. It will then be made explicit. 

A new generator of evolution @ determines the time dependence of the 
elements of the vacuum while the correlated elements are determined by a 
creation operator C acting on the elements of the vacuum. The computa- 
tional rules of these operators are well known and only the relevant part for 
our purpose wil be recalled later on. 

Since the emitted photons cannot interact with the electron (by construc- 
tion), the kinetic equation for fss'sS'oo"' '^^^^ involve only elements 

fss'si'nn'' ^^'^^ (with possible other values for the indices) where 

the r) indices arc of the incident i type. We label by an index inc such clc- 

ments and since by definition, the field described in f f ^ 

has never been in interaction with the Dirac particles, we can factorize the 
particle and field dependence in these functions: 

~^^l...^ls^l'^,■■■n[^^l...ils^^'g,■..^i'l , . _ rAti---MsAi^,...M'iMi---/isAj,---M'i /,w /,x o\ 

Jss'ss'nn' \'') ~ Jss'ss'OO {'l^)Jnn'inc\.'^) {°-^) 

fnn'inc{t) describes the incident field and can be of quasi classical or arbi- 
trary, as considered in previous sections. Its time evolution is independent 
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of the evolution of the Dirac field and corresponds to the propagation of a 
free field. Of course, the resulting field depends, via the constitutive relation 
( |8.2D , also of the emitted field and of their interference. Keeping that factor- 
ization in mind, the kinetic equations describes now a closed equation 
for the Dirac correlation functions, in canonical momentum variables. The 



relation ( 5.10 ) and the similar ones for more than one particle allows to 
transform these kinetic equations into kinetic equations for the reference 
distribution functions, with the help of relations ( |6. 13 -6. 16). 



From Q and ( ^TTl) , we know the link between the Wigner distribu- 



tion functions (in variables position and mechanical momentum) in terms 
of the correlation functions in the other set of particles variables and the 
electromagnetic field. Through the creation operator, the correlation func- 
tions involving the self-field are determined by the vacuum elements only. 
Using the factorisation (^.3|), the Wigner distribution functions can finally 
be expressed in terms of the reference distribution functions. Therefore, the 
closed evolution equation in one set of variables (the reference functions) 
can be transferred into a closed evolution equation for the Wigner functions 
in the other set of variables. If the gamma's ( |5.9D are known, the kinetic 
equation for the Wigner functions is completely determined (and is depen- 
dent of the free incident field that has to be considered as given) and it is 
possible to check if divergences are present. 

8.2 Explicit link between the mechanical and canonical mo- 
mentum distributions 

We make explicit the link introduced in the previous subsection between 
the Wigner distribution function and the correlation functions for the Dirac 
particles. 

From ( |5.9D , the doubly Fourier transform one article Wigner distribution 
fW( 



functions f^{—K., —X) and the correlation functions are connected by: 



1^1, ■■■'^n u'...u', 



/„^(-K, -X) =h^T.J d^'^i • • • / d^K J d'k[ ... 1 d^k'^, Yl E 

Y.jd'p j ^V7^^'^'"'^'(P,p',fci,...fcn;fei...fe;o/ri 

^^^J.' 



llOOnn' 
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+ / '^V / (P'' P ' ^1' • • • ^1 • • • K')foi01nn' 



The 7's are defined through (|5.9D in Lorentz gauge. The /'s functions are 



defined in ( 4.1^ ) with a convention that the arguments are those of ( |8.1| 
when they are not exphcitely written. 

We use the constitutive relations (18.21) to obtain for instance: 



/noorm'(P' P'; kii/i; k2i^2; • • • ; k„i/„; k'^i^^; ^1/3; . . . ; k^'f^/; 

= H /noonn'(P' P'; kii^iryi; k2i^2f?2; • • • ; KfnVm 

iv} 

k;i/;r?i;k'2i/^ry^;...;k'„,i/;,r/;,;f) (8.5) 

Considering the invariance of the distribution and of the relation (^]^) with 
respect to the permutation of photons, in each term of the sum Yl'{n}' 
can place all photons of the same type (i) in the last indices, with a counting 
factor. Obvious conventions are in use for values of m and m' near their 
limits. 

/noOnn'(P' P'; kii/i; k2Z^2; • • • ; k„i/n; k'lZVi; k'^u!^; . . . ; k^/t-^,; t) 

= E 



ml{n — m)\ m'\{n' — m')\ 

^/noonn'(P' P'; kii/io; . . . km,i/„o; km+ii^m+ii • • • ; k„i/„i; 

;t) (8.6) 

Considering the equivalence conditions, the distribution functions with an 
index o are numerically equal with the distribution functions where the index 
o is replaced by s. That operation has to be considered as realized when 



( |8.5[) is reported in (8^). We have then the factorization similar to 
between the incident field and the electron with its dressing photons. The 
index i is useless for the function fnn'inc where all photons are of that type. 

/noonn'(P' P'; kiJ^i; k2Z/2; • • • ; k„zy„; k'ii/[ ; k2i/2; . . . ; k^,^^'; 



E 



I m\{n — m)l m'\{n' — m')! 
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></nOOmm'(P' P'; kil/is; . . . kml^rnS; k[u[s; . . . Ka'^'m'S] t) 

^ f(n-m){n'-m')inc(km+ll^m+l ■ ■ ■ ^knyn'-i^rn' +l^rn' +1 ■ ■ ■ jkn''^n'!^) 

(8.7) 

In the subdynamics, the correlated components (involving e.g. the self-field) 
are obtained by the action of a creation operator C on the elements of the 
vacuum and we have: 

/noomm'(P' P'; kii^is; . . . k^z^ms; k'^i^^s; . . . \i!^,v'^,s] t) 

= ^ < limmm' ^i^JL'\C\ss's~s' > fss'ss' (8.8) 

ss'ss' 

where the spin indices have not been written. 

Dynamical restrictions are present in the summation J2ss'ss'' the parity 
of the sum s + s' + s + s' and the value of 5 = s — s' — s + s' remain invariant 
under the evolution. We place accordingly a prime on the summation symbol 
to remember that property. 

Therefore, combining the two previous relations, we obtain: 

/iToonn'(P' P'; kiz^i; k2Z^2; • • • ; k„ ; k'^ i/[ ; ^1^2; • • • ; K''^'n'',i) 

n! n'l ' 

= ~~u — '' TT ,u 7TT H < 1100mm'^/i'|C|ss'ss' > fss'Ss'oo 

^-^ , m\[n — m)\ rn'\[n — m )\ 

m,m' V / \ / ss'ss' 

^ f{n-ni){n'-m')inc0^m+l^m+l ■ ■ ■ '^^nl^n^^m' +l^'m' +1 ■ ■ ■ ^n' '^'n' ^) (8-9) 

Inserting that relation into ( ^.4| ) leads to: 

/^(-K, -X) = ^3 ^ 5^ 1 d'h . . . / d^K I d'k[ ... 1 d'k'^, 

ss'ss' n,n' 

~~n — ' TT /,/ / " ^ < lWOmmn^j,'\C\ss'ss > fss'sS'oo 

, m\(n — ml! m'\(n — m'y. 

m,m' \ / \ / 

^ f{n-m){n' -m')inci^m+l'^m+l ■ ■ ■ ] ^n'^n] ^m'+l^'m'+l • • • ' ^'n'^'n''^^) 

fljl' 
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~~w — ' TT /,/ / " TTT < 001lmm'ilfL'\C\ss'ss' > fssTs'oo 

, m\[n — my. m'lln' — m')! 

m,m' \ / \ I 

^ ~F7 — ^ ^ ,u ,' TT, < OWlmm'n'fi'\C\ss'ss' > fss'ss'oo 

^ , m\{n — ml! m \[n — m l! 

m,m' \ / \ I 

^ f{n-m){n' —'m')inc0^m+ll^m+l • • • ! j^iU^^i . . . ; k^/Z/^/; t) 

Ef}} 7l!\ ~ 
~~n TT /,/ / " TTy < 1010mmV/i|C|ss'ss' > /ss'SS'OO 
, mlln — ml! m'\(n' — m'y. 

m,m' \ / \ I 

^ f{n-m){n' —m')incO^'m+ll^m+l ■ ■ ■ ] ^n'^^n] ^m'+l^m'+l • • • ) ^n'^n''i ^) 



(8.10) 



8.3 The dressing operator 

Similar relations can be written for all the joint Wigner distribution func- 
tions involving an arbitrary number of particles. They can be symbolized 
with the use of a dressing operator x~^{finc) Hi; depending on the incident 
field distribution functions, as (r is the number of particles described in /^) 
(the definition is general although we use it for one particle only): 

f^= J2 < r\x~Hfinc)\ss'ss' > L'ssW (8.11) 
ss'ss' 

The elements of x~^(/mc) for r = 1 can be read directly on expression ( ^.10| ) 
for the choice of variables used in that expression. Since the correlation func- 
tions that do not involve the field can be replaced in terms of the reference 
function (see section 6), the dressing operator could also be defined as the 
link between the reference and Wigner functions. 

That relation can be inverted, for the relevant values of the indices 

noted by a R, to provide: 

fssTs'ooR = E < ss'ss'R\x\r > (8.12) 
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From now on, we will leave out the writing of that index R: a restriction on 
the relevant values of ss'ss' is implicit. 

The dressing operators X considered above present new features 

with respect to the dressing operators considered up to now (Q and refer- 
ences within) . First of all, it is a time dependent operator since the incident 
field distribution fine evolves freely with time. This point should not appear 
as unexpected. Indeed, the canonical and mechanical momentums of the 
electrons-positrons are linked in a naive approach by transverse components 
of the potential vector A that evolves with time and contains obviously also 
the external field and not only the self-field. 

A second difference is that, in no coupling limit (q^ — > 0) of the dressing 
operator, such as it is defined in ( |8.11 ) and ( 8.12| ) with respect to the corre- 



lation functins, the dressing operator does not reduce to the identity. If we 
define it with respect to the reference function, it reduces to the identity in 
the no coupling limit. 

In Q, in a quantum optics context, we have constructed the dressing 
operator by imposing a form of the kinetic operator so that some observ- 
ables are naturally defined in the physical representation, among them, the 
observable linked with the probability of finding the atom in the upper state. 
Here, the dressing is necessary to deal with the physical variables in place of 
the auxiliary ones (canonical momentums) that are required for the quan- 
tization of the Dirac field. The two points of view are closely connected. 
In the present approach, we may begin with an a priori proposal for the 
gamma's (|5.9| ) or obtain conditions on them by imposing the finiteness of 
the evolution operator. 



8.4 The evolution equation of the Wigner distribution func- 
tions 

We know that fss'sS'ooit) satisfies a close kinetic equation in the subdynam- 
ics: 

dtfss'ss'ooit) = ^ < ss'ss'\e{finc)\rr'rr' > frr'ff'oo{t) (8-13) 

where the operator depends on time via the incident field. The dependence 
of the operator in the spin variables is not written explicitly. From the 



37 



connection ( ^.11 ), we have successively 

dtf^it) = E < r\dtx-\nncmss'~s~s' > L'ss'ooit) 

ss'ss' 

I 

+ E < Ax'^{finc{t))\ss'ss' > dtfss'ss'OO 
ss'ss' 
I 

= ^ <r\dtX'^{finc{t))\ss'sS > f ss'ss' OQ{t) 
ss'ss' 
I 

+ E < Ax'\fincmss'~S~s' > 
ss'ss' 
I 

X ^ < ss' Ss'\Q{finc{t))\rr' ff' > fj-r' ff'ooit) 
I 

< r\dtX-\Lcmssss' X ss'ss'\xiLcit))W > /f W 



« ss'ss' 
I 



+ E E <r\x'\hncmss's§' > J2 < ss's~s'\e{nncmrr'fr' > 

X < rr'fr'\x{fincmu > fu it) (8.14) 
Therefore, 

dtf^it) = E < r\mnc{t))\u > f^{t) (8.15) 

u 

with 

< rmnnc{t))W > = < r\ {dtX~Hhnc{t))) x(.Lc{t))W > 

+ E < r\x-Hhnc{t))e{f,nc{t))x{finc{t))W > (8.16) 
u 

That set of equations (8.15) is the basis of the dynamical equations of QED 
in the physical variables.. 

As will be established soon, the kinetic operator Q{finc{t)) presents the 
usual ultraviolet logarithmic divergence. Our aim in this paper is to exam- 
ine whether $(/jnc(^)) presents the same illness. We will not proceed to a 
systematic analysis of that kinetic operator $(/inc(i)) but we focus on the 
contributions that are divergent in 6(/j„c(^)). They appear already at the 
second order in the charge via the self interaction. Although our approach 
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is based on an expansion in powers of the charge, we do not reahy assume 
the analycity of the operator. We are aware that resummations of class 
of contributions may be required, as well known in statistical mechanics in 
similar context ||8|, ^ for dealing with non-analycity. For the first orders we 
treat here, no problem should arise. 

For examining the divergence property, we may restrict ourselves to con- 
tributions that do not involve the incident field. Indeed, the incident field 
provides a natural cut-off at high energy in its wave number expansion. 
Therefore, we consider the system of equations ( ^.15 ) for initial conditions 



where no incident field is present. The influence of the self- field on the 
spreading of the one particle distribution function can still be studied under 
these initial conditions. Such kind of problem is not usually considered in 
QED while it is treated as a first application of the Schrodinger equation in 
non relativistic theory. 

8.5 The zero*'* and first order evolution operator 

The zero*'' order < 1| |1 > operator is provided by: 

< >=< l|x-^(°)e(0)x^°^|l > (8.17) 

The free motion operator 0*^"-* is identified with Cd given in (6.3). It does 
not involve transitions in the numbers ss'ss of the distribution function and 
is a diagonal operator in the canonical momentum representation . 

The zero*'* order of the operator < l|(x~^)'^'^^|ss'ss'00 > can be read in 
( |8.1[1| ) by considering the values n = n' = m = m' = {). The relation will 
be inverted, by computing the inverse of a 16 x 16 matrix, to provide the 
operator < ss'ss'00|x^'^^|l > (see later on ( ^.28|) ). 

The first order < 1|<I>(^)|1 > operator involves the incident field and is not 
considered here. Such contributions are present for describing the absorption 
of the incident free field but the contributions of self-energy retained here 
requires at least a the second order in the charge. To evaluate the order of 
the contributions in equations ( 8.16| ), we have to take into account that the 



creation operators C in ( |8.1ClD are at least of order m + m' in the charge. 

Since for instance 7^/3^™' involves a factor q^~^^ , the first correction to 

is of the second order and arises from the contributions where m = 1, m' = 



or m = 0, m' = 1. (x "^)*-^'' is thus obtained from ( ^.10 ) by considering these 



two contributions while taking the creation operator at the lowest order. 
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8.6 The second order evolution operator 

The second order < 1|<I>(^)|1 > is thus provided by: 

< 1|$(2)|1 > = < l|(;^-l)(2)eWx(°)|l > + < l|(x-^)(°)e(2);^(0)|l > 

+ <1|(x-^)Wg(0)x^')|1> (8.18) 

We determine first of all x*-^"* in terms of (x^^)^^^- The relation XX^^ = ^ 
provides the condition X^^Hx^^)^'^^ + X^'^Hx^^)^'^^ = 0; hence we have 

The kinetic operator ^^"^^ can thus be written as: 

(8.20) 

In order to analyse the possible singularities of we can consider the 
operator t^(o)$(2)(-y~1)(o) for which we have: 

;,(0)${2)(^-l)(0) ^ 0(2) + [;,(0)(^-l)(2),e(0)] (8.21) 

The operator S^^) is an operator in variables position and canonical mo- 
mentum. From the subdynamics, we know that (P and Q = 1 — P are the 
projectors respectively on the states of the vacuum and of the correlations, 
the usual L is ihC, with a similar link for PQ^^^P with the usual collision op- 
erator 

^^(2) g (PG(2)p = due to the lack of factors in the definition 

(iH)) 

ihPe^^'^P = ihPC'QC^^'^P (8.22) 

The explicit computation of a relevant element of 0*^2) be found in 
appendix D. It will be further analysed later on. It contains the usual 
logarithmic divergence for its part diagonal in variables K and P. That 
divergence has to be compensated by the commutator in the second term 
of ( 8.21| ) to provide a finite result. That compensation can be studied as 
well on ( |8.18| ) as on ( |8.21| ). The advantage of ( 8.21| ) is that we work in 
the original variables and we can study separately the consequences of the 
dressing. 
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8.7 Expliciting the zero*'* order dressing operator (x ^)*^°^ 

We focus now on the second term in ( ^.211) and determine first of all the 
(^-i)(0)^ ^(0) operators. If we call /J(°)(-K, -X)(= /i°)(-K,-X)) the 
links hand side of ( |8.10| ) when no incident field is present and the {x~^) 
operator is replaced by {x~^)^^\ we have: 



/„7)(-K,-x) = n^Y. J d^p J rfV7^^;^^'°°(P,P')/rfo'ooo 

W T.Jd'pJ d'pjlf^'\p, p)/ro^ooo (8.23) 

f^g^^ (— K, —X) can be identified with the reference function /^^^ (— K, — X)). 
We take the value of the relevant 7's in the expressions ( 6.13 ), ( |6.14D , ( |6.15|) 
and ( |6.16| ). In order to obtain an algebraic relation in place of ( ^.231) , we 
replace the 7's by their expression and take the Fourier transform with re- 
spect to the variable X. All the dependence in X lies in exponential and the 
Fourier transform provides Dirac delta functions. We multiply both sides 

by ^^^p e R and integrate over X to obtain: 

(-K, P) = E / dh J dy 5{v - p' + hK)5{Y - P - ^K) 



xe;''^(p',P)/rfoooo(P,P') 



E / J d^P Hp -p' + hK)6{P + p' - ^K) 



x^P'(-P,-p')/&o(P,p') 



I d^P I d^P Hp + P - KK)5{V + p' 



-K) 

2 ^ 



xCp' (P,-p')/oVioo(p',P') 
+^^E / J d^pS{p + P + hK)6{P - p - ^K) 
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x5^^'^-p,p)/roU(P,P) (8-24) 
Performing the trivial integrations, we obtain 

/„T(K,P)= E €f^(P-^K,P + jK)/W(K,P) (8.25) 
where the matrix d^p is given in (|6l^ ) and the non vanishing elements of 



the one particle distribution function /^^^'^'(K,P) are by definition: 
/W(K,P) = ft3/|^/^J^^^(p + ^K,P - ^K) 

/W(K,P) = -r.V£oo(-P + ^K, -P - ^K) 
/W(K, P) = h'f^g,,,{V + ^K, -P + ^K) 
/W(K,P) = fi3/^;^^oo(P - ^K, -P - (8.26) 



From ( |8.25| ), /^^"'^'^'(K, P) is indeed associated with '^^ and ( ^.24 ) implies 



the association with /foiooo- Similarly, /^^^/^' (K, P) is associated with (J^'^ ^'^ 



and (8.24) implies the association with /qiqioo- 

This new set of functions {/«'^«'^'(K, P)} can replace everywhere the set 
of functions present in the right hand side of the previous expression and the 
elements of the dressing operator can be defined as the connection between 
/^(K,P) and /«'^«'/^'(K, P). The relation (|]2|) allows to identify in these 



variables the new elements of the {x '^)^^^ dressing operator with those of 
the 5 matrix: 

{X-')T'''' = - P + ^K) (8.27) 

8.8 Computing the zero*'^ order dressing operator 

The dressing operator {x~^)^^\ linking /^^°^(K,P) and /?'^?'^''(K, P) (and 
thus also the four functions in the right hand side of ( 8.26| ) , is thus diagonal 



in the variables K, P and this property has to be shared by x^^^ ■ We can 
determine the form of x*^*^^- Using the expression of 5, the unitarity of r, the 
invariance of the trace, the explicit form of our basis matrices B^p as direct 
product and the properties of the Pauli matrices, x(°) is easily seen to be: 



X':r = -Af (P + P - 2 K) (8.28) 
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8.9 Free motion operator 0*^°^ 

Let us provide the one particle free motion operator Q^^^ ( ^.17 ), (7.3) in 
the representation provided by ( 8.26| ) for the Wigner distribution function 
/^(K,P). The diagonal operator Q^'^^ has the following non vanishing 
elements: 



P-|K 



< l^l^'KP|G(°)|l/il/i'KP >= — (£^p+iLK - Ep_nK 

ZJt, V 2 2 

< i^VKP|e(°)|i/ii^i'KP >= ^ (^_p+|K - 

< lulu'KP|e(°)|l/il/i'KP >= — (e^.hk + E_-p,!tK 

< l^l^'KP|e(°)|l/il^i'KP >= — (-Ep_^j^ - E_p_nj^ 

ill, V 2 2 



i.29) 



8.10 Second order diagonal motion operator 6^^^ 



From the appendix D, we have for the four contributions corresponding to 
the four diagrams (the lower indices refer to the vertices of appendix C 
involved) : 



< KPI 



KP >=i cPk 



X < cos 



2 |P+ir»K| 2 \'P+\r{K~h\L\ 



2eo(2^)3wk 



■ cos 



+ sin 



\P+khK\ 



■ sm 



2 ^|P+i;iK| ^|P+i;iK-;ikl 
-3 cos sin 3 cos 



2 2 
^|P+iRK-Rk| . , ^\V+lhK\ 



sm 



+ 



P + IkK -Kk P + IkK 



\P + \KK - fik\' \P + \hJ^\ 



sm I 



|P+^nK| ^"^'^IP+ifiK-fikl 



~-^|P+i;iK-;ik| + -^|P+iriK| ~ 
For the second diagram, we have: 



(8.30) 



< KPI 



7.1 



KP >= i / d^k 



2eo(27r)3wk 



X < cos 



"T COS ^- h sin - 



iriK-;ik| _ |P-i?iK| 



■ sm 
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2 |P-i;iK| |P-l;iK-;ikl \P~^hK^hk\ |P-lnK| 

—3 cos = sin = 3 cos = sin = 

2 2 2 2 

/ P - IhK -hk P + IhK \ 1 

+ T • T sin6',„ 1 sinT 



1 



-^|P-iRK-Rk| ~ -^IP-ir^KI + 
For the third diagram, we have: 

< KP|<^^|KP >= J d'k^^ 



(8.31) 



2eo(27r)3a;k 



an Of) 
2 |P+1rK| 2 \P+htiK+hk\ 2 IP+i^^K+Rkl 2 |P+^?^K| 

X < 3 cos cos h 3 sin sin 

2 2 2 2 

2 ^\P+lhK\ 2 ^\P+hhK+h-k\ 2 ^IP+^'^K+Rkl 2 ^|P+i'^K| 

— cos = sin = cos = sin = 

2 2 2 2 

/ P + iftK + hk P + \ . ^ . ^ 1 

V|P + ^ftK + nk| |P + iftK|; \P+^m |P+^;.K+;.k| J 

^ (8.32) 



-^|P+inK+nk| + -^IP+iftKI + 
For the fourth diagram, we have: 



< KP|0i^,i^,p7^|KP >= -i f d^k — 

I A/|a3f|a25|i J 2eo(27r)3a;k 

2 V-^?.K| 2 V-^'^K+nkl V-lnK+nk| . 2 V-^nK| 

X < 3 cos ^ cos h 3 sm sm 

2 2 2 2 

2 ^|P-inK| ^IP-inK+nkl ^IP-inK+nkl ^|P-lftK| 

cos = sin = cos = sin = 

2 2 2 2 

P - \KK + h'k P - \KK 



V|P-^ftK + nk| |P-^ftK|; \P-^m |P-^;.K+;.k| 

1 



P-iftK+ftk| ~ -^|P-iriK| ~ 



(8.33) 



From ([4.12| ), we have sin^p = Let us consider first the problem of 

convergence of the integral over k. The problem of convergence could arise 
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for large wave numbers k. In that hmit, we may use the properties: 

hm sm^ = hm cos"' = - (8.34) 

fc— *oo 2 fc— >oo 2 2 



arising from 



2 9p Ep + mc^ 1 mc^ 

cos = = o + (8.35) 

2 2i?p 2 2ii'p 

. 2 -Ep - mc^ 1 mc^ 

si'^ o ~ „ = - — 777;- (8.36) 

2 2Ep 2 2E/P 



We look in appendix E how the various integrands behave in that limit. 
Combining partial results, we obtain 



I{'^k)M\pr + li^kjhAp + I{^k)M\Tr 

q^c^ m?c^ I 1 



127) 



2eo(27r)3a;k2fi2^2 ^E^^^^,^^^ ^|P-inK| 
For P and K in the non relativistic domain, we obtain 

I{^k)M\prNR + H^k)M\FNR + H^k)M\TrNR 

g^c^ 1 -fiK.P 

~ ~^ 2eJ^(27r)3^2fi2^ m 

= -i^^(-)3-^c^ / ^ (8.38) 

and the mass correction 6m is 

5m= I d^kS^c^ } „ =3^ ^ r dk^^ (8.39) 

J 47reo 87r2?ia;^ 47reo 27r /ic Jo fe'^ 



We examine whether the second term in the right hand side of ( p.21| ) 
provides similar divergences for the proposal (^.8|). Details of the compu- 
tation are provided in the appendixes. The result can be synthesized as 
follows: in the non relativistic domain, small K and P, we obtain obviously 
a behaviour in K.P with a coefficient that diverges also logarithmically with 
a cut-off for large wave numbers: the same kind of behaviour as for the col- 
lision operator B^^^. Therefore, we can ascertain that the transition to a 
description in terms of the mechanical momentum provides new terms that 
are comparable to the divergent part for the description in terms of the 
canonical momentum. 
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9 Conclusions 



This paper treats the interacting Dirac and electromagnetic fields as a statis- 
tical system in the sense of a A^-body system, N arbitrary large. Therefore, 
we have introduced an analogue of the reduced distribution functions for 
the determination of the mean values of the observables. Our main preoc- 
cupation has been the obtention of the distribution functions describing 
the variables position and mechanical momentum associated with the elec- 
tron. These functions can be interpreted as the physical representation |Q] 
describing the observables and these functions have to replace the original 
description in terms of the canonical momentum. The transition between 
the two sets of functions is ensured by a dressing operator previously intro- 
duced in the context of quantum optics. 

A main requirement on the dressing operator is that it leads for the phys- 
ical representation to evolution equations free from the usual divergences. 
The usual trick to eliminate of these divergences is through counterterms to 
provide mass (charge, . . . ) renormalization, assuming an infinite bare mass 
for the charged particle. However, we think that the origin of the divergences 
lies in the use of unphysical variables (canonical momentums) that have to 
be considered only as useful intermediate entities. Progress along these lines 
requires progress in the association of operators describing physical observ- 
ables with classical functions of position and mechanical momentum. The 
present paper is a prerequisite to tackle those questions still pending. The 
two aspects (choice of the association and finiteness) can enlight each other. 
No contradiction with the renormalization program is expected: the latter 
may serve as a useful guide for the determination of the dressing operator. 

The main differences of the present work with usual QED lie in inter- 
pretating differently the origin of the divergences and in providing a more 
satisfactory conceptual framework. We hope that it will induce further work 
on the subject. 

A The charge and current operators 

The explicit computation of the charge and current operators requires the 
evaluation of the product of the spinors Up_^ and Up_^ by the adjoints. For 
instance, the coefficient of c^(p)c^/(p') in the expression of pn is given by: 

^*P,,iUp',f^' = uo^,f*{p)T{p')uof,> = (r+(p)T(p'))i^i^, (A.l) 
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The explicit evaluation of all terms is lengthy but straightforward. In terms 
of elements of the Pauli matrices and of the angle i^p,p' ^^^^ satisfies the 
fundamental trigonometric relation on the sphere: 



Gp Opi . 6p . Op' ^ ^, 
cos cos + sm sm -^p.p , 



cos (yCp,p: 
we get for the coefficients a: 

«slt;.'(P>k) = (cOS<^p_l^k,p+ifikV 



(A.2) 



/((pxk) 



XT,,,,' sm ■ 



sm ■ 



//(pxk) 



-^hk\\p + ^nk.\ 



. ^Ip+lfikl . f^lp-iftkl 

.an' a sm sm ^ 



(A.3) 



(A.4) 



p + \Kk 



p + Uk| 



2-- 1 ''lP-5^k| . f^lp+ifikl 
' cos sm — 



- ' — ) cos — ^ — si"^ — ^ 



^ip+i?ik| . f^ip-ifiki 



|p- i?ik| 



,11 

S flfl 



(P,k) 



-a 



.,ii_ _ 

S fj,fj, 



(P,k) 



(A.5) 
(A.6) 



11 / 1 \ le*^ 



P + ^hk _ ^|p-iftk| 
|p + ifiki 



cos ■ 



sm ■ 



p - ^Hk 



sm ■ 



^|p-iRk| ^p+ink' 



p - ifik| 2 
/ / p + i;ik 



cos 



MM 



lp+ifik| 



xcr^ ^/ cos ■ 



'|p-ifik| . ^|p+iftk| 
— sm 



p — ^hk 



sin ■ 



COS ■ 



(A.7) 



11 / 1 \ 9eC 



|p + ifik| 



COS ■ 



sm ■ 



?ik 



'^hk\ 



sm ■ 



COS ■ 



'MM 
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+ 




COS ■ 



sm ■ 



|p+irik| 



sin ■ 



cos ■ 



IP+i'ikI 



(A.8) 



11 / ^ \ ^(^^ 



a^f,, cos V2p-i;ik,p+iRk 



P 



|p + i;ik| |p-ink| 
p + ifik 



sm ■ 



^|p-irik| . ^|p+iftk| 



P 



|P - 



|p + iftk| 



MM 



X sm ■ 



'|p-iftk| 



■ sm ■ 



Ip+i^ikl 



■ sm ■ 



p + ^hk \ p 



IP ■ 



(A.9) 



a 



'(P,k) 



11 _ 

/ik 



SkA- 



-il, 



cosv9p_i^l,,p+iRk 
p \ . ^|p-irik| . ^|p+iftk| 



tj.fl' 



|p + i;ik| |p - ifikl 



p + iftk 



sm ■ 



■ sm ■ 



'MM • 



P 



^hk\J |p + iftk| 



p + irik \ p 



iftk 



MM 



|p + irik|^ 



IP ■ 



iftk| 



. '^|p-i;^k| . f^lp+iRkl 

X sm sm 



"A«M'(P'l^) 



aAuM'(P'l^) 



(A.IO) 



(A.ll) 



The elements of a with a subscript minus can be btained by replacing in 
the above expressions e\^)^ by e.j^^^ 



B Elements of the joint characteristic operator 

The computation of the elements 7^^'' ^ "'^'^(p, p') and 7^^^ ^ ''^(P; p') is straight- 
forward. We need the first contribution in the expansion of ( |6.1| ) in powers 
of qe- We have obviously: 

g-H.X.A,,(r)^X.V+«.r„ ^-H.X.A,.{r) 
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X.V+iK.r 



B. 



(B.l) 



Introducing the relation ( p. 21) , we obtain 



g-H.X.AxKr)eX.V+iK.r^^^g-H.X.Ax.(r) 



:-^X.A^Kr)e^^K.X^.K.rgX.V5^^ 



(1) 



(B.2) 
(B.3) 



i/ii^t 10 (|p^ p'~) jg given by 



i/^Vio 



ip'.r 



d're J* uo/.''r(-p') 



xe 2 



iK.Xg^K.rgX.V 1 gj^^^^,Xe*l-ri?,^e^r(p)noM 



(B.4) 



Acting with the displacement operators e-^'^, we obtain: 

1 ' 



|K.X^.K.r_^^^^^^_Xe^k,.(r+X)^^^^ 



xe 2^ 



ip.(r+X) 



T"(p)^iOM (B.5) 



The integration over r can be performed and provides a Dirac delta function. 



7„/3 (P>P)- 



a/3 



(5(p-p' + K + ki)e^^K-^ 



(B.6) 



7^i;^^'°^(p,p')--^ 



1 



h'^ h ]j 2eo(27r)3 
:ek;vrXe --^ '^(p,p) 



<5(p-p' + ftK-Sk;)e 



/^„^iK.X 



(B.7) 
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lap (P'P) 



1 —iqe / h 

2eo(27r)3 

»p'.X 



(5(p - p' + /iK + ;iki)e^2 



iK.X 



x-^ek,.,Xe^^f^'(-p,-pO 



(B.8) 



lap (P'P) 



1 — ZQe / ^ 



5(p - p' + ?iK - ^ik'Je 



/^„^iK.X 



X ^ek;., .Xe-'^i-^e^ 5^^^' (-p, -p') 



(B.9) 



7^^/'^'^°(p,pO = 



(B.IO) 



7r'^''^^(p,pO 
1 



1 —iqe / /' 

2£o(27r)3 



6(-p - p + /iK - /ik;)e*2 



/AjiK.X 



(B.ll) 



7a^ (P,P) 
1 



1 —iqe / ft 
F n Y2eo(27r) 
X 



r5(p + p' + ;iK + ftki)e^2 



iK.X 



ek...X(e^'^-^-l)e^if'^(-p',p) 



(B.12) 



7S'"'°'(P'P') = 



(B.13) 



C The evolution equations for the correlation func- 
tions 



This appendix provides all the non-vanishing contributions due to the part 
CiD of the evolution operator for the correlation functions. Let us first 
consider the contributions due to the part jCjdi of Cjd involving the product 
of creation and destruction operators c^ca. The contributions diagonal in 
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the occupation numbers ss'ss' can be written as: 
< ss'ss'nn'\jCim\ss'ss' nn' — 1 > 



M i=l,s' j=l,n' 



^ 2£o(27r)3u;k, 



2 ^' 



xs{p'i + hk'j,ii;p[,fi',)l n 5(k;+i,z/;+i;k;,z/;) 

yr=j,n'-l y 

Two contributions correspond to the same matrix element of JCjdi- 
< ss'ss'nn'\CiDi\ss'ss'n + In' >„ 



(C.l) 



xa 



11 



(Pi + -k„+i)S'(p^ + ?lk„+i, /x; p^, /x-) (C.2) 



< ss'ss'nn'\C[Di\ss'ss' n + In' >b 

l^n+l Lt' 1=1, S V ^ 



X -/ii/.'(Pi - -fik„+i,-k„+i)5(pi - hkn+i,n';pi,fii) (C.3) 

The contribution involving a change in the occupation numbers ss'ss' can 
be written as: 

< ss'ss'nn'\CiDi\s + 1 s' + Iss' nn' — 1 > 



E 



h 



-,^2£o(27r)3a;^. 
J d^pal}_^^,{p, -'k'j)S{p - ^hk'j,ii';ps+i,iis+i) 



;5(p + ^?ik;.,M;p;'+i,/x',/+i) ( n 5(k;+i,z^;+i;k;,z/;) |(c.4) 

\r=j,n'-l 



We now consider the contributions due to the part jCim of jCjd involving 
the product of creation and destruction operators c+ca"^. The contributions 
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diagonal in the occupation numbers ss'ss' can be written as: 
< ss'ss'nn'\CiD2\ss'ss' n — In' > 



fx' i=l,s j=l,n 



2eo(27r)3u;k, 



xS{Tpi + h'kj,ll']Y>i,IJ'i)\ n 'S'(kr+l,fr+i;k^,fr 

\r-=j,n-l 

Two contributions correspond to the same matrix element of Ci d2 ■ 

< ss'ss'nn'\CiD2\ss'ss'nn' + I >a 



(C.5) 



E ^K'-.i)E/^^C-,i E . 2.0(2.)'.,, 

M »=1,S \ ^ n'+l 



11 



,(pi + ^fik;,,+i,k;j,+i)S'(pj + ?ik;,,+i,/i';pj,/Xi) (C.6) 



< ss'ss'nn'\CiD2\ss'ss'nn' + 1 >b 

The contribution involving a change in the occupation numbers ss'ss' can 
be written as: 



< ss'ss'nn'\C[D2\s + 1 s' + 1 ss' n — 1 n' > 

\s+s'+s+s' 



1 



^H'Y^i-ir^'^'^^' y si. 



MM 



j=l,n 



2£o(27r)3a;k, 



ir=j,n-l 



We now consider the contributions due to the part Cms of Cm involving 
the product of creation and destruction operators 6+6a. The contributions 
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diagonal in the occupation numbers ss'ss' can be written as: 
< ss'ss'nn'\CiD3\ss'ss' nn' — 1 > 



1 

ih 



^^EE E 



M i=l,s'j=l,n' ^ ^ J 



x5(p: + ;ik;,/x;p:,/i^ n .5(k;+i,z.;+i;k;,z.;) (c.9) 

\r=j,n'-l I 

Two contributions correspond to the same matrix element of Lioz'- 
< ss'ss'nn'\CiD3\ss'ss'n + In' >„ 

= EE/ ^^^n.i E t/ 2.o(2vr)3., 

I^n+l M 1=1, s' V ^ ^ «-n+l 

x«i!+i-M/i^(Pi + J/^k„+i,-k„+i)S(p^ + /ik„+i,^;p^,/x9 (C.IO) 



< ss'ss'nn' |>C/£)3 1 ss'ss' n + In' 



n+l 



X"i.„+i-/ii/.'(Pi - -^k„+i,-kn+i)S'(pi - hK+i,n';Pi,fii) (C.ll) 

The contribution involving a change in the occupation numbers ss'ss' can 
be written as: 

< ss'ss'nn'lil/DsIss' s + ls' + lnn' — 1> 



MM 



xS'(p + ^/ik^-,/x;p5,+i,/i'g/+i) I Jl 5(k(.+i,z/^+i;k^,f^; 

kr=i,n'-l 



(C.12) 



We now consider the contributions due to the part CjDi of Cm involving 
the product of creation and destruction operators 6+6a+. The contributions 
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diagonal in the occupation numbers ss'ss' can be written as: 
< ss'ss'nn'\CiD4\ss'ss' n — In' > 



ix! i=\,s3=\,n V ^ "3 



\r-=j,n-l 

Two contributions correspond to the same matrix element of Lib^'- 
< ss'ss'nn'\CiDi\ss'ss'nn' + 1 >a 



(C.13) 



ih V n+i^z^y 2£o(27r)3(jk' 



1 



n' + l '^"^ ^ 



< ss^ss'nii!\CiD4\ss'ss'n'n! + 1 

= i^-^EK-;'+i)E/^X'-ME, 



n'+l 



\ 2eo(27r)3a;k/ 

^ n'+l 



xa 



11 



;(Pi - 2^^n'+i>k;,/+i)5(p^ -;ik;,,+i,M;p^,Ai) (C.15) 



The contribution involving a change in the occupation numbers ss'ss' can 
be written as: 



< ss' ss' nn'\CiDA\ss' s + Is' + In — In' > 
1 ' 
ih 



X 



2eo(27r)3a;k, 

j c?^paJ.,VM,.'(P,kj)S'(p + ^/ikj,/i';p5+i,/i5+i) 



^•^Ip - ;^fikj,/^;p's'+i,/i5'+i) I n 'S'(k^+i,t'r+i;kr,t'^ 

ir=i,n-l 



(C.16) 



We now consider the contributions due to the part Cid^ of Cid involving 
the product of creation and destruction operators c^h^a. No contribution 
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is diagonal in the occupation numbers ss'ss'. The term corresponding to 
the diminution of two occupation numbers is: 

< ss'ss'nn'\CjD5\s — 1 s' s — 1 s' n + In' > 



i—l+s—j 



^n+l i=l,S j=l,S 

1 

2' 



X I n 'S'(Pr-+l,/ir-+i;Pr-,Mr 

vr=i,s—l 

xS{-{pi + pj),f„+i;kn+i,z/„+i) 



yr=j,s-l 



(C.17) 



The terms corresponding to a transfert of occupation numbers are: 
< ss'ss'nn'\Cjr)5\s — 1 s' s s' + In -\- In' > 

-l+s+s' / 



l^n+l n' 1=1, S 



2£o(27r)3a;k, 



xa^;^^^_^.^/(pi - -?ik„+i,-k„+i) ]J S'(pr+l,Air+i;Pr-,Air) 



I r=i,s— 1 



(C.18) 



< ss' ss' nn'\Ci Db\s s' + Is — Is' n + In' > 



2£o(27r)3a;k„, 



+1 



11 



\ r=i,s— 1 



xS'(-pi + ftk„+i,/v,;p',,+i,/^^,_^J (C.19) 
The term corresponding to the majoration of two occupation numbers is: 
< ss'ss'nn'\j0.jD5\s s' + Is .s' + 1 11 n' — 1 > 



2£o(27r)3a;k, 
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X J d3paiJ_^^,(p,_k^)5(-p + ink^.,/x';p^-,+i,/ig'+i) 



x5(p + ^?ik^-,M;Ps'+i,/^','+i) I n '5(k;+i,f;+i;k;,i/;; 

Lr=j,n'-1 



(C.20) 



We now consider the contributions due to the part jCid6 of jOid involving 
the product of creation operators c^b'^a'^. No contribution is diagonal in 
the occupation numbers ss'ss' . The terms corresponding to the diminution 
of two occupation numbers are: 

< ss'ss'nn'\jC,iD6\s — 1 s' s — Is' n — In' > 

= -- E E E <-.)(-ir^+^'-^ ' ^ " 



1 



\r=i,s—\ 

XI Jl 5'(pf+i,/xf+i;pf,/if:) Jl S'(kr+i,fr+i;kr,t'r 

>r=j,s-l / \r=;,n-l 



(C.21) 



< ss'ss'nn'I^C/Del^ — Is's — Is'nn'+l > 
1 

ih 



~ E E E (-1)-^^^-. 

i=l,s j=l,s 



n' + l 



^ 2.o(27r)3a;^(p^^^^.) 



Xa^}^^^+^.-.(-(pi-Pj),--(pi+Pj-)) I n '5(pr+l,Mr-+i;Pr,Mr-) 



ir=i,s— 1 



\f=j,s-l I 

{C.22) 



The terms corresponding to a transfert of occupation numbers are: 
< ss'ss'nn'\CiD6\s — 1 s' s s' + In — In' > 



■l;.3^.(.0EE(-ir^^'^'' 



l=l,n fj,' i=l,s 



2£o(27r)3wk, 
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yr=i,s—l 

\r=l,n—l 



(C.23) 



< ss'ss'nn'\jCiD6\s — 1 s' s s' + Inn' + 1 > 



ir=i,s— 1 



xS(-pi - /lk;,,+i,/i';p'5,+i,/i'5,+i) 



(C.24) 



< ss'ss'nn'\CiD6\s s' + Is — 1 s' n — In' > 



-l+s'+s' 



l=l,n M i=l,s 



2£o(27r)3a;k, 



xS{-pi- hki,n;p'g,^i,lj,'^,^^) { Yl S{'kr+i,i^r+i;K,J^r 

\r=l,n—l 



(C.25) 



< ss'ss'nn'\Ci£)Q\s s' + 1 s — 1 s nn' + 1 > 



^ 2eo(27r) 



n'+l 



xa 



11 



,i(-Pi - 2^'^™'+l'^n'+l) I n •S'(Pr+l,/ir+i;Pr, Ar 

\r=i,s—l 

xS{-pi - h\c'^,_^i,tJ,;p'^,+i,fi'^,^i) (C.26) 
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The term corresponding to the majoration of two occupation numbers is: 
< ss'ss'nn'\CiD6\s s' + 1 s s' + In — In' > 



MM 



j=l,n 



2£o(27r)3cc;k. 



\.r=j,n-l 



(C.27) 



We now consider the contributions due to the part Cjdj of Cjd involving 
the product of destruction operators bca. No contribution is diagonal in the 
occupation numbers ss'ss'. The terms corresponding to the diminution of 
two occupation numbers are: 

< ss'ss'nn'\CiD7\s s' — lss' — lnn' — l> 



^ E E E (-1)-^^^-'- 

i=l,s' j=l,s' l=l,n' 



x'5(k| + i(p^ + p;.)W n '5(p;+i,/i;+i;p;,/i;)J 

\T=i,S —1 J 



\f=j,s'-l 



1 n S{K+iA+iX,K)\ (C.28) 

r=Z,n— 1 



< ss'ss'nn'\CiD7\ss' — lss' — ln + ln'> 



i-l+s'-j 



^n+l i=l,s' j=l,s' 



\ 2eo(27r)3u;^(p,+p,) 



xa 



ll+i-fi'.n'.iliPi-'P'j)^liPi+'Pj)) [ n '5(Pr+l,^r+i;Pr,/^r^ 



Lr=i,s'— 1 



XI Yl '5(Pf+l,/if+i;Pf,/if) 'S'(--(Pi + Pj),t'n+i;k„+l>i^n+l) 

^f=i,s'-i / 

(C.29) 
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The terms corresponding to a transfert of occupation numbers are: 
< ss'ss'nn'\CiD7\s s'— Is + ls'nn'— 1> 

l=l,n' n' 1=1, s' V ' 



Xal}_^,^,XPi + -hk'i,-k'i) ^ n 'S'(Pr+l,Mr+i;Pr,Mr; 

^r=i,s'— 1 



xS'(-p^ - ;ikJ,/x';Ps+i,/is+i) I n 'S'(kJ.+i,t'^+i;k;,f^; 

kr=/,n'— 1 



< ss'ss'nn'\jCiD7\s s' — ls + ls'n+ln'> 



(C.30) 



n+l 



yr=i,s'— 1 

xS'(-p- - ;ik„+i,ju';ps+i,/xs+i) (C.31) 



< ss'ss'nn'|>C7£)7|s + 1 s' s s' — Inn' — 1 > 
1 



i=l,n' M i=l,s' V ' '^i 

\r=i,s'— 1 



I r=(,n'— 1 



(C.32) 



< ss'ss'nn'\CiD7\s + ls'ss' — ln+ln'> 
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\r=i,s'—l 

xS{-p[- hkn+i,n;ps+i,fis+i) (C.33) 
The term corresponding to the majoration of two occupation numbers is: 
< ss'ss'nn'\CiD7\s + ls's + Is' nn' — 1 > 



X / c^^pai|-^^'(P,-kj)S'(-p- -/ik^-,)U;Ps+i,As+i) 



^r=j,n'-l 



(C.34) 



We now consider the contributions due to the part Cjds of Cjd involving 
the product of creation and destruction operators 6ca+. No contribution is 
diagonal in the occupation numbers ss'ss'. The term corresponding to the 
diminution of two occupation numbers is: 



< ss'ss'nn'\CiD8\s s' — 1 s s' — 1 nn + 1 > 
1 

ih 



i=l,s' j=l,s' 



\ 2^o(27r)3a;i(p,+^,) 



„ii 



n 'S'(Pr+l' /^r+l! Pr' /^r, 
\r=i,s'—\ 

X -^Cl (Pi + Pj- ) > i^n'+l ; k^'+i ^v'n'^l) 



(C.35) 



The terms corresponding to a transfert of occupation numbers are: 
< ss'ss'nn'\CiD9,\s s' — 1 s + 1 s' nn' + 1 > 
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xa 



11 . ' 1 



Pi ~ 9^k^'+l' '^n'+l) I H 'S'(Pr+l) A*r+i; Pr) /^r. 



yr=j,s' — 1 

x5(-p- + ;ik^,^_i,/i;ps+i,/is+i) (C.36) 
< ss'ss'nn'|£/£)8|s + 1 s' s s' — Inn' + 1 > 



xa 



11 



1 



'^'(-Pi + 7J^k^'+li k^'+l) H 'S'(Pr+l,/ir+i;Pri/^r) 



\ r=j,s' — 1 

X'S'(-p- + ftk^/+i,;u';ps+i,;Us+i) (C.37) 
The term corresponding to the majoration of two occupation numbers is: 
< ss'ss'nn'\CiDs\s + 1 s' s + 1 s' n — In' > 



X y (i^|3ai^^+^^,(p,kj)S'(-p + ^ftkj,/i;p5+i,/i5+i) 

xS{p + ^h'kj,fi';ps+i,fJ's+i) i n 5'(kr+i,M,.+i;k^,t'r) |(C.38) 

D The diagonal second order kinetic operator 
e(2) for one electron 

The operator 0^^^ is the operator in variables position and canonical mo- 



mentum ( |8.21 ). We will be interested in its elements that connect state 
1 110000 > to itself, with the same value for the arguments of the functions. 
These terms allow to illustrate obviously the problems associated with the 
self-energy interactions. The basis for their computation is the relation 
( |8.22| ). The terms of interest to us are: 

< 11000019(2) IllOOOO >= — < 110000|«?iP£'QC(^) 1 110000 > (D.l) 

The explicit forms of the last vertices are found in appendix C and have to 
be combined with the expression of the C^^^ operator. As for all elements 
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of the subdynamics operators, we begin from: 



27rz Jc z — iC 



(D.2) 



and we have to compute the residue arising from the vacuum states only to 
obtain 

Since at first order, 



1 



(1) 



1 



-iC- 



1 



^z — iC J z — iCo z — i£o' 
the first order creation operator has a very simple form: 



(D.3) 



< 110010/ii/i'i|C(^)|110000/i2/U2 > 

_ -1 

~ 27ri 



/ dz < 110010/zi/i'J —iCiD2 — 1110000^2/^2 > 

Jc z — iCo z — iCo 



(D.4) 



where the prime on the integral sign means that only the contribution of the 
pole corresponding to a vacuum state has to be included in the integration 
path (here, only the second propagator satisfies the criterion). Only jCid2 
provides a non-vanishing contribution. 

Let us list all the required matrix elements: The combination of vertices 
which provide an ip diagonal in the variables are: ( p.3 /C.5), ( |C.7| / p.lD , 
( C.24 / C.32 ), (IC.31 / C.25 ). To each combination is associated in the usual 
way a diagram If we explicit the variables (in the K, P representation), 
we need the following expressions, written in Lorentz gauge, of appendix D 

takes four values: three values for A and one value s). From ( p.3| ), ( p.5| ), 
O, O, don, (^, (^, (^, we obtain 



< P + ^hK, fi;P- ^hK, n'; 0; 0; 0; 0\Cidi 

1 1 

|P + -hK - hk, fl;P- -hK, /x'; 0; 0; k, i^; > 



In' ^ 

ih Y 2eo(27r)3wk 



a 



.11 

U — flfl 



1 



1 



(P + -hK - -hk, -k) (D.5) 



< P + ^hK - hk, fl;P- ^hK, /x'; 0; 0; k, u; 0\£id2 
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1 1 
|P + -/iK,/xi;P- -hK, fi' ■,0;0;0;0 > 

< P + ^fiK, /x; P - ^/iK, m'; 0; 0; 0; Q\Cid2 
|P + ^/iK, P - ]^KK - hk, fi; 0; 0; 0; k, >b 

< P + ^/iK, P - ^hK - hk, fl; 0; 0; 0; k, i^\Cidi 
\P + ^hK,n;P- ^?lK,//'i;0;0;0;0 > 

= — — a}} „',-7(P- ^^K- ^;ik,-k) (D.8) 

< P + ^hK, /i; P - ^/iK, /x'; 0; 0; 0; 0\jCid6 

|0; P - ^hK, n'; 0; -P - ^HK - hk, fi; 0;k,u > 

< 0; P - ^hK, /i'; 0; -P - ^hK - hk, fi; 0; k, u\£id7 
|P + ^;iK, /ii; P - ^/iK, /i'; 0; 0; 0; > 

< P + ^/iK, /x; P - ^fiK, m'; 0; 0; 0; 0\£id7 

|P + ^hK, n; 0; -P + ]^KK - hk, fi; 0; k,u;0> 
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< P + ^hK, /x; 0; -P + ^hK - hk, fi; 0; k, u; 0\Cid6 
|P + ^hK, /i; P - hK, 0; 0; 0; > 



1 

ih 



h 11 l.„ 1. 



. X, -a,i„/,.(P- -?iK + -ftk,kXD.12) 



Equations (5.39) provide the relations between the a's and the (5's, for / s 
The denominators which appear in the four diagrams are, after the in- 
tegration over z has been performed: 



1 



1 



ie + E,-p_i 



-^IP-ifiKI + ^'^k 



P-iriK-rik| ~ -^|P-ifiK 
1 



"a.24/C.32 



-^C.3l|/]a25 



ie + -Eipii 



^IP+ifiKI + ^^k 



P+iRK+Rk| ^ ^|P+iRK| 
1 



ie - -E^|p_ifiK+Rk| ~ -^|P-ifiK| ~ ^^k 



(D.13) 
(D.14) 
r (D.15) 
(D.16) 



In all propagators, the ie's play no role since the denominators never vanish. 

Inserting the value of the different elements, we obtain for the expression 
of the relevant elements of C*-^-*: 



< llOOlO/ii/i'ilC^i) 1 110000^2^2 > 

= Hn' 
-I r' 1 



^ 11 , IfcU u \xKr 



2eo(2vr)^"^^+'^^^^^P+2'''^''^^)'^' 



dz ■ 



1 



^(p + nki;p) 



h 11 / 'cKr 



-1 



dz 



2TTi Jc Z- i^p + lEp. -UJk.Z- i£;p+;,ki + i^p' 



■5(p + riki;p) 



2eo(27r)3cuk, 
1 



'^Ai+fii^2 



(p + ^ftki,ki),5, 



M'i,M2 



■S{p + ftki;p) 



(D.17) 
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It can easily be realized that -Ep+^jki ~ ~ ^"-^fci cannot vanish, so that no 
prescription involving a ie is required. Similarly, 

< 110001^i/i'i|C(^)|110000/i2^2 > 

_ s3 / ^ ^,11 , IfciJ ^J\xKr 



5(p' + ;ik;;p') (D.18) 



< 010101^1 /i'l I C(^) I IIOOOO//2/U2 > 



ft3 / ^11 (_p'_l;jk'i,-k'i)(5fV 

^ 5(-p'-ftk;;p) (D.19) 



< 101010/ii I C(^) 1110000^2/^2 > 



r—T^ oF-u,,' „' (-P - -^ki, ki)(^^ 



1 



S{-p-hki;p') (D.20) 



Let us now focus on the spin dependence of the diagrams which is pro- 
vided by the product of the vertices only: the propagators play no role in 
that respect. We take into account the expression of the a's for the computa- 
tion of the spin dependence of the diagrams and we introduce new functions 
S"s. Their definition takes into account the compensation of the factor ^ 
in ( p.E| - p.l"2 ) by the factor ih present in ( |D.1| ). In a similar way, the factor 



in ( p.5| - p.l"^ is compensated by the factor in the a's. The factors 
arising from the a's are taken into account in the expressions of the 5"s. We 
keep explicit the s^v) factor. 



flu 

V 
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The coefficient of the spin identity operator will be called 
of the Pauli matrix S^^^^. In a similar way, we have: 



and that 



(D.22) 



For the last two diagrams (|C.24| / 1C.32| ) and (|C.31| / 1C.25| ), it ought a priori 
to be taken into account that the label of the argument present in the spin 
operator (present in the expression of the a's) requires a spin component 
opposite to the index of a. Since we have to perform an summation over 
the spin variables, this point does not play a role and the contribution of 
the two last diagrams is similar to that of the first two diagrams. 

We will see that only the scalar parts 5^/^, 5^^^, S^^ /c^, 
"^csi/^^i '-'^ "^'^ vanish. From the expression oi uie coetticients 

a, we obtain by direct identification, taking into account the arguments re- 
quired in ( D.5 - |D.12|) : (we shall consider successively the expression with the 
scalar and with the vectorial photon.) 



a 



11 

S fj,fj,' 



'(P,k) 



cosV'p-ink.p+ir.kV 



h{pxk) 



p — iftk||p + iftk 



■.<r^^/ sm- 



■ sm 



a 



.11 

S flfj.' 



/ , N QeC 



a 



11 

S flfl' 



(P,k) 




|P + ^ftk| 



cos 



COS ■ 



|p+i?ik| 



■ sm 




/I /I' 



(P,k) 



(D.23) 
(D.24) 



with the convention fj,' = —fi'. From those expressions, we deduce, taking 
into account the value of the arguments of the a's in ( D.S - DTT^ ): 



-QeC cos fp+^fiK-h'k,P+^hK 



ft((P + ^ftK)x(-k)) 
P + \KK - hkWP + ^KK 



sm ■ 



\P+\hK-hk\ 



sm ■ 



(D.25) 
2 

(D.26) 
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(D.27) 
(D.28) 



The contribution to the scalar part 5 j^r^^ becomes therefore: 



qy cos2 93p+ir,K-r.k,P+ir.K 



+qtc h # :Tr; tt^ — sm sm 

^ |P+ ifiK- ftk|2|P+ i;iK|2 2 2 

(D.29) 

while the contribution to the vectorial part clearly vanishes considering 
(P.27I - P.28D . 

The corresponding contribution to the scalar part 
tained by changing the sign of K: 



7^ can be ob- 



(«0sags + b^-bg;!!,) = 9eC cos V'p- ir.K-ftk,P-isK 

2 2,2 fc'|P-^^Kp-((P-iftK).k 
|P-UK-?ik|2|P-UK|2 



2„2.2fc'|P-^^K|2-((P-iftK).k)2 . ^V-ifiK-nkl . 2V-i^K| 



■ sm 



■ sm 



2 

(D.30) 



while the contribution to the vectorial part clearly vanishes. 

The computation of the contributions due to the combination of vertices 
( 1021 / 10321) and (|03l| / p35| ) requires the expressions: 



(D.31) 



and: 



b |a24| s = geC 



P + i/iK + ;ik 



|P + i/iK + ^ik 



cos ■ 



nP+iRK| 



■ sm ■ 



^|P+i;iK+;ik| 



P + iftK 
|P + iftK| 



'^|P+i;iK+;ik| . '^|P+l;iK| 
cos = sm 



C.32S 



bc.24s 



(D.32) 
(D.33) 



>27s = <leC 
P 



P - iftK + nk 
|P - \KK + nk| 



"^IP-iriKI f^lP-iftK+Rkl 
cos ^ sm 



f^lP-inK+nkl "\P-^hK.\ 
cos sm 



(D.34) 
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C.25s 



(D.35) 



The contribution to the scalar part of S^T^IFT^ and 5^r^xpr^ requires the 
computation of: 



'.243 



32s 



2 2/ 2 '^|P+5'iK| . 2 

c cos sni 



IP+isK+nkI 



+ cos 



2 ^|P+i?iK+Rk| . 2 ^|P+5?iK| iF + lhK + hk P + iftK 

~ 2 \P + ^HK + hk\' \P + lhK\ 



sm 



2 2 
X sin^ip^i^KI sm0|p+i;,K+;.k|) 



(D.36) 



Ic.sik-tfesk = I cos 



2„2 



2 ^|P-iftK| . 2 ^|P-i?iK+Rk| 



sm 



+ cos 



|P-isK+nk| 



sm 



2 2 
2-|P-inK| iP-iftK + nk P-iftK 



2 2 
X sin ^ip. 1 ;jK| sin ^|p_ i ;,K+;.k|) 



2 |P - iftK + riknP + iftK| 



(D.37) 



Anew, the contribution to the vectorial parts of 5 vanishes. 

We now consider the contributions involving the vectorial photon. It is 
useful to write the expressions of the a in the following way, in order to 
allow a more direct identification of the a's and the b's: 



;p,k) 



-SkA- 



P-iftk . „, 
+ - f- rr sm cos ■ 



^|P-iRk 



P + ifik ^|P-i;ik| '^|P+i;ik| 
■ cos sm 



\P + ^hk\ 



iftk| 



.QeC 



- ekAX- 



ekA^^ 



2 
P 



MM 



;ik 



cos ■ 



|P-irik| 



■ sm ■ 



f^lP-inkl f^lP+irikl 
sin ^ — cos 



|P+irik| 



(D.38) 



Tik 



|P + i;ik| |P-iftk| 
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^|P-iftk| ^|P+iftk| 
X sin sin ■ 



2 2 

ekA cos V^p-ifik.P+irik 



ekA- 



|P-Mk| / |P + i?ik| 



P + iftk / P + Uk 
+ hkA-T; ' 



"^IP-irikl '^|P+irik| 
X sm sm 



(D.39) 



"A+M/.'(P''^) = -^T^^MM'^kA- 



X sm sm ■ 



?ik 



X- 



^^^^^ 



2 2 

ekA cosv9p_i^j^_p^i^^j^ 



ekA- 



P - iftk \ P + i?ik 



iftkl / |P + iftk| 



+ ekA-- 



P + \h\i \ P 



p + i^ikiy |p - iftk| 



"^IP-irikl "^IP+irikl 
X sm sm 



(D.40) 



with the relation ( |A.2| ) 

From those expressions, we derive: 



aC.3A = leCG-^X- 



|P + i/iK| 



COS ■ 



IP+inK-nkl 



+ 



P + iftK - nk ^|P+ifiK-;ik| 
|P + i?iK - nk| 



sm ■ 



cos 



■ sm - 



^P+|kK 



|P+iftK| 



ac.3A 



(D.41) 
(D.42) 



and therefore, the corresponding contribution to the first term of 
can be obtained as 



a^a||x = ^eC^ cos 



2 P+iriK 



sm 



2 ''|P+|?iK-Rk| 



+ sin 



X sin 6* 



P+|rK _2 ^|P+f iK-^ikl 1 P + \hY^ - nk P + i^K 
2 2|P+ i/iK-nkl'lP + ^nKl 



cos 



(D.43) 
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The computation of the contribution of the second term of >-'p^/ 
the expressions: 



,5 requires 



= iqeC 



P + iftK - hk 
P + iftK - hk\ 



^|P+iriK| ^|P+ifiK-?ik| 
cos = sm 



ekA^ 



P + iftK 

\p + ^nK\ 



sm cos = 



>2feA 



Using the relation: 
we obtain: 

A 

/ P + IhK 



Yl (ekA^c) . (ek;,xc') = 2c.c' 
A 



(D.44) 
(D.45) 

(D.46) 



P + ^hK - hk 
IP + IhK - hk\ 



IP+inKI |P+i;iK-fik| 
cos sin 



ihKl 



sm ■ 



2 

IP+^'^KI 



cos 



|P+^riK-rik| 



(D.47) 



The expressions ( D.43| ) and ( p.47 ) determine the contribution due the vecto- 
rial photon to the scalar part 5|^riyPr^ while the contribution to the vectorial 
part vanishes considering ( p.42D -( p.45D . 

We now consider the contribution of the vectorial virtual photon to the 
second diagram ( p.l /C.7). From 

„ fP-^hK ^|P-inK-nk| ^|P-i;^K| 
= Qecei^x- \ — f-^cos ^ sm- 



JP + i;iK| 

P-^hK-hk ^|P-inK-Rk| ^P-^nK 

+ - f sm cos 

|P - i/iK -hk\ 2 2 ^ 



(D.48) 



o^A = (D-49) 
the corresponding contribution to the first term of 5^/5^ can be obtained 



as 



7^0^ cos^ ^ sm 



P-inK . 2 ^|P-isK-?ik| 
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+ sin 



X sm ( 



2 P-I'-'K 



COS 



2 ^|P-iriK-rik| _^ 1 P - - hk P - ^^K 



2 IP 



j/iK-nkl IP 



|P-^riK| |P-^riK-rik| 



itiKl 



(D.50) 



In a similar way, the computation of the contribution to the second part of 



C.l/C.l 



requues: 



cos = sm 



P - i/iK 
|P - ^hK\ 



sm ■ 



|P-iRK| 



■ COS ■ 



|P-ifiK-rik| 



from which we derive: 



P - ^hK 
IP - iftKl 




cos ■ 



sm 



sm ■ 



|P-^nK| |P-iriK-rik| 



o 2 2 I 2 2'^i^\ . 2 

zg^c cos sm 



cos 



(D.51) 
(D.52) 



|P-^ftK-fik| 



ihK\ 



+ sin 



cos 



P-^ftK-ftk 



fftK -hk P 



|P - ^hK - hk 
The contribution to SI? 



sm ( 



IP- 



fRK-nkI 



(D.53) 



c 7|/|cj| ^^^^ obtained from ( p. 50]) and (p.53 ). It can 
be obtained also from the expressions for the diagram ( |C.3| / |C^ ) by changing 
the sign of K. Anew, the contribution to the vectorial part vanishes. 
We now consider the contribution of the vectorial virtual photon to the 



third diagram ( p.24 /C.32). The computation is similar and we obtain: 

f hk 



0'C.24\ 



X sm ■ 



■ sm ■ 



P + iftK + nk| \P + lhK\ 
^|P+iriK+rik| 



(D.54) 



71 



from which we deduce: 



o loii A = -Q|a24|A (D.55) 
A:2|P+ i;iK|2-|k.(P + i;iK)|2 



2 2fc2 



|P + iftK + nk|2|P + l?iK|2 



X sm 



2 



sin — 



2 2 
The second part of 5^"^lFr^ requires the evaluation of 



b|a2|x = QeC [ekACOS(^p+i;jK,P+iaK+rik 
P + ifiK \ P + ^ftK + 
■|P + ifiK| ] |P + iftK + nk| 



P + i/iK + Tik \ P + iftK 



P + ^hK + hk\l |P + iftK| 



. "^IP+i/iKI . 
X sm sm ■ 



'|P+inK+nk| 



(D.56) 



(D.57) 



br 



C.35A — QeC 



e^A COS V'p+iftK.P+iyiK+^k 



ekA 



+ ekA-- 



P+ iftK 
|P + iftK| 



P + iftK + hk 



P + iSK + hk 



P + iftK + hk\ 



X sm ■ 



■ sm ■ 



|P+ ^hK + hk\ 
P + iftK 
■|P + iftK| 

'|P+iriK+rik| 



(D.58) 



and therefore, we obtain: 



E b|al]Ab|32lx = 3qy cos^ 1 ,k,p+ 1 ,K+.k 



-4ggC^ cos 



(P + IhK + hk) P + i/iK 



P+^nK,P+^nK+nk \^ |p + + nk| ■ |P + ^hK\ 



X sin = sin ■ 



fTiK+nkl 
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ekA- 



P + \kk 



P + \KK + Kk 



P + \KK\ / |P + ^KK + Kk\ 
P + \tiK + fik \ P + ^^K 



X sin 



P + \fiK + Kk] / |P + 2^K| 
P+^riK+?!,k| 



|P+^riK 



sm 



o 2 2 2 
O^gC COS if. 



P+i;iK,P-i;iK+;ik 



-4ggC^ cos ip. 



/ (P + iriK + nk) P + ifiK 
P+^nK,P+lnK+nk l^|p+ Ir^K + nkl 'iP + ^/iKl 



^|P+ir.K| 
X sin = sin ■ 



fftK+ftk| 




P + \riK P + iftK + Kk 



P + ^?iK|'|P + ^ftK + nk|_ 



X sm 



|P+^riK| 
= sin 



|P+^riK+?ik| 



(D.59) 



The contribution due to the vectorial photon to 5^r^lJ^r^ is obtained from 
( p.56| ) and ( p.59 ) while the contribution to the vectorial part vanishes. 

We now consider the contribution of the vectorial virtual photon to the 
fourth diagram ( p.31 / Cl2^ ). From: 



hk 



P - \KK 



X sm ■ 



|P-inK| 



sm ■ 



-\KK + hk\ |P 

-iriK+rik| 



«a.25A 



we obtain: 



(D.60) 
(D.61) 



2 2fc2 



PIP 



Ik.fP 



iftK)|2 



X sm 



= sm 



P-inK + ftk|2|P-i?iK|2 
^|P-iRK+Rk| 



(D.62) 



73 



From: 



br 



as^x = qec ekACOS(/Jp„i^K,P-ir>K+rik 



ekA- 




^hK + hk\ 



X sm ■ 



P - ; 

|P-|fiK 



HK + nk\ / |P - ^hK\ 



sm ■ 



(D.63) 



C.25A 



ekA COS '^p_inK,P-ir»K+ftk 



ekA- 




P - i/iK 



ftk 



+ ekA- 



P - iftK + hk 



X sm ■ 



P - i 

2 



HK + nk| / |P - UK| 







sm ■ 



|P-i?iK+rik| 



we obtain: 

5Ib|a3^-b^2|A+ = Sg^c^cos^ V^P-i^K.P-iftK+r.k 



(D.64) 



4qlc'^ cos (/'P-i^iK 



2,.^^,P-iriK+rik 



'(P-i^K + nk) P 



;ik| IP 



X sm ■ 



|P-inK| 



;riK+/ik| 



2 

1 + 



sm ■ 



P - i?iK P - i/iK + ?ik 



\P - ^hK\' \P - ^hK + hk\ 



X sm 



■ sm 



[P-inK+Rk| 



(D.65) 



is obtained from 



2 2 

The contribution due to the vectorial photon to ^^T^TT-^gs 
(A.a75) and (A.a78) while the contribution to the vectorial part vanishes. 

We now collect our partial result to obtain the complete contribution for 
each diagram. For the first one, we have to add the contributions of ( p. 29] ), 
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( p.43| ) and ( p.47| ) and we obtain, for the scalar part S^T^j^ ^ of Sc.-j^f^ 
(A.a29), taking into account the value of s{u): 



-qlc^ cos^ V'p+ir.K-rilcP+iriK 



2 2,2 fc'|P + I^Kp-((P + ifeK).k)2 . 2 

— (/gC ^ — ^ r-^^ ,„,_ — ^-i — TTTT^ — sm 



IP+iaK-nkI 



2 2 



COS 



P + i/iK- ;ik|2|P + i;iK|2 
hK.\ , 2 "|P+^?»K-?ik| 



sm 



2' 

P+1 



■ sm 



+ cos 



-^nK-Rkl 



^2 ^| P+|?»K| 
2 



sm 



riKl 



+ 



1 P + ^ftK V + \KK 



2 IP + IftKl • IP + IhK - hk\ ^"^'^IP+l'^K-.kl ^in'^IP+i.KI y 

/] 

|P+iriK-rik| 



2„2 



cos 



|P+^riK| 
= sin 



^- h COS 



■ sm 



|P+^nK| 



1 P + ^hK P + ^hK 



hk 



sm ( 



2|P + iftK|"|P + i;iK-nk| 

2 2 2 

- -q^c COS V'p+irjK-rik.P+lriK 

2 fc^|P + iftKp-((P + ifeK).k)2 

"'^'^"'^ |P + i;iK-;.k|2|p + i;,K|2 



|P+i;iK-;ik| ^"^'^|P+i;iK 



IP+iriK-nkI 



sm 



2 ^|P+iaK| 



2 2 
-9eC 



3 cos 



|P+^rK| 



sm 



3 cos' 



2 2 

^|P+iriK-rik| . „ ' 



■ sm 



1 P + iftK P + ^hK 



hk 



sm ( 



|P+i;iK-;ik| ^"^'^|P+ifi,K| 



2|P + iftKnP + i;iK-nk| 
Using the explicit form of cos(/5p p;, we obtain: 

V+^'^K-Rkl 



a|/|ai 



2„2 



cos 



2 ^|P+i^ 



COS 



|P+^r»K| 
2 

(D.66) 



+ 



P + iftK 



nk P 



hK 



1 P + iftK - hk 



P + iftK 



P + iftK-nkl'lP + i/iK 

1 

2 

IP - 2'-l - ^ 2 



sm 



'|P+iriK-rik| 



■ sm 



|P+ir»K| 



2 IP + iftK - nk| • IP + i?.K| '^"^|P+I''K| sinf.|p+i,K^nk| 



2 2t2 



fc2|P + ifeKp - ((P + ^ftK).k) 2 2 ^|P+ir»K-ftk| 2 ^|P+|'iK| 
|P+ iftK- ftk|2|P+ iftK|2 



■ sm 



■ sm 
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2 2 I o 2 |P+|^K| . 2 V+lnK-nk| 2 IP+^'-^K-r.kl . 2 |PH 

-gg c I o cos sm ^ h 6 cos ^4 sm 



1 P+|ftK P + ^hK 



sm6'|p+i^K-ftk| sint^|p+i;iK| 



Let us note that we can write: 



|P + i;iK + ftk|2|P + inK|2 
|P + i^K + hk\^\P + iftKp - |(P + iftK + hk).{P + iftK)|2 



i;iK + rik|2|p + i;iK|2 



1 



|(P + inK + nk).(P + iftK)|2 
|P + i;iK + fik|2|P + i;iK|2 

so that we obtain: 



(D.67) 



5r 



2„2 



COS 



hK\ 2 ^|P+|?»K-Rk| 



■ cos 



(D.68) 



+ sm ^ sm 



\P+^hK\ 



, 2 2 
+9eC 



2 2 y 

2 ^iP+inKi V+^y^K-Rki 

3 cos = sin = 



+ 3 cos^ 



|P+iriK-rik| 



■ sm 



|P+^rK| 



P + ^hK P + iftK 



hk 



\P + ^hK\' \P + ^hK - hk\ 



smf 



|P+iriK-rik| ^"^^^IP+ifiKI 



(D.69) 



We now consider the contribution of the second diagram which can be ob- 
tained by changing the sign of K in the previous expression: 



C.7/C.1 



2 2/ 2 ''IP-^'^KI 2 ^\P~^r^K^r^k\ 

7„ c cos cos 



. 2 ^\P-lnK-nk\ . 2 
+ sm sm ■ 



^2 2|o 2 |P-^^K| . 2 V-^nK-nk| 2 |P-^^^K-r.k| . 2 V-^^K| 
+gg c I 3 cos sm ^^^^ h 3 cos ^4 sm 



i/iK P - iftK - nk 



^;iK| |P - ^hK - nk\ 



sm^^|P_i,^K-rik| sint^|p_iftK| 



(D.70) 
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We now consider the contribution of the third diagram which can be ob- 
tained from (pIS^ ), ( pT56|) and (pl59|) : 



2 2/ 2 

C COS 



^|P+ifi.K| . 2 ^|P+i?iK+nk| 



■ sm 



+ COS 



2 V+llK+ftkl . 2 



sm 



1 P + iftK + P + i/iK 



2 IP + IftK + hk\ • IP + i;.K| ^"^'^IP+I'^KI ^"^'^IP+l.K+.kl 



I 2 2t2 



fc^lP + I^Kp- |k.(P + ^ftK)|^ 
|P + i;iK + ?ik|2|P + iftK|2 



„ \P+lhK\ ^\P+lhK+hk\ 
X sin = sin = 



-Aqlc^ COS 1 



(P + ^hK + hk) P + ^hK 



X sm 



P+^nK,P+^nK+nk \^ |p ^ + nk| ■ |P + ^hK\ 

a Q 
|P+^riK| . |P+^RK+Rk| 



2 

1 + 



■ sm ■ 



P + \riK P + \KK + Kk 



|P + ^riK| |P + iftK + nk|_ 



2 |P+iriK| '^|P+ifiK+fik| 
X sin = sin = 



(D.71) 



Using the exphcit form of cos and the relation ( D.6^) , we obtain 



■^ 10.241 / 10.321 — ~9eC^ ( ^OS 



■ sm 



+ cos 



2 ^\P+\hV.+h-k\ . 2 ^\P+^hK 



sm 



S^gC^ I cos 



2 f^lP+ir^KI 2 '^IP+l'iK+ftkl 



■ cos 



2 ^|P+iriK+rik| o IP+^iK] 
- sm sm 
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2 2 P + i^K + ftk P + iKK 



iP + ifiK + ftkl'lP + i^Kl 
X sin ^|p+i;,K| sin ^|p+ i r.K+r.k| 

We now consider the contribution of the fourth diagram which can be ob- 
tained by changing the sign of K in the previous expression: 

■^ 10.311 / 10251 - ~^eC- I COS ^ Sm ^ 



2 '^|P-i;iK+;ik| o '^|P-i;iK| 
- COS sm 



2 



22/ 2 ^\p~\m 2'iPzi^K+^ 

-SggC I cos ^ cos ^ 

^|P-iriK+rik| o ^IP-ifiKI \ 
h sm"' ^ ^ sm.'^ 



„ „ P - i/iK + h\<i P - i;iK 



|P - i;iK + ^ikl'lP - \kk\ 
X sin0|p_i,^K| sme|p„i,,K+r.k| P-^S) 

We know therefore all the required elements of the operator 0^^). The 
presence of a saclar part only for describing the spin dependence of the 
computed elements of that operator shows that they are also diagonal in 
the spin indices, and not only in their variables K and P. 



E The logarithmic divergence of the collision op- 
erator 0^^) 



The dominant part of the integrand in these expressions are (for the part 
that does not vanish through angular integration): 



1 



1 



dom 2eo(2vr)^a;jj Ihuj^ 



1 



dom 2eo{2'Kyuj\^2hoJk 



(E.l) 

(E.2) 
(E.3) 
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-^(^fc)M |a3l| / |c.25| 



1 



(E.4) 



dom 2eo(2vr)^C(;k 2hujk 

The sum of these contributions clearly vanishes. We now consider the cor- 
rection due to the propagator. 



-^('^fc)M |a3l / |a5 



Pr 



En 



-k. 



|P + iftK| 



(-l)c(P + -K).ksin0|P+i,K| 



(E.5) 



-^(^fc)M|a7|/|a 
P 



+k. 



-T P 



h 



2-^2 (^|P-i/.K| 



-K).ksin0,p_i,^K| 



(E.6) 



-^(^fc)M |a24| / |c.32 
P+ i/lK 



Pr 2eo(27r)3u;k4?i2a;2 



2-2 (-^|P+i?iK| 



+k. 



\P + ^hK\ 



c(P + -K).(-k)sin^|P+i,K| 



(E.7) 



g2c2 1 / 

= -^2.o(2^)3^k4fi2.;2 i-^|P-i.K| 
+fe- |plf|;K| ^-lMP - ^K).(-k) sin V-i.K|) 



g2c2 1 / 

(io^ 2^ Wp+^'-^KI - ^|P-i.K| 



(E.8) 



Therefore, 

, (k.(P + iftK))2 . ^ (k.(P - i;.K))2 . ^ ^ 

+c — ; r : — sm fip , i ^k-i — c — ; r : — sm t^ip i sx^-i 

|P + iftK| IP+a'iKI |P-iftK| |P-2?^K| 



(E.9) 



We now consider the correction due to the k dependence in the trigonometric 
functions. We use ( |8.35| ) and ( |8.36| ).The last terms, in sin 0jp_|_i^j^|do not 
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contribute. 



mc 



Tr 2eo(27r)3u;jj 



X 4 cos' 



4 sin' 



^\P+^hK\\ 1 



-2hu!k 



(E.IO) 



MCl/CA 



die' 



mc 



Tr 2eo(27r)3a;k 2fta;k 



X 4 cos 



|P-inK| 



4sin^ 



(E.ll) 



X 4 COS 



.32 



Tr 
inK| 



mc 



\P+^hK\ 



2hiOi^ 



(E.12) 



-^(^A:)Af |a3l| / |c.25 



ay 



mc 



Tr 2eo(27r)3u;j^ 2?icjk 



X 4 COS 



2 ^|P-|/^K| 



4 sin 



2^|P-i;.K|\ 1 



-2huj}^ 



(E.13) 



We consider the corrections due to the /c-dependent factors. We have 



(P + hiK - h\^f = h'^k'^ - 2(P + ^hK).hk + ... 



1 



1. 



h'k' ( 1 - 2— (P + -?iK).k 



(E.14) 



Therefore, 



|P + -ftK -hk\ = hk(l 



1 



1 



\P + ^hK-hk\ hk 



i(P + l«K).k)+.„ 



(E.15) 
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Using these relations, we have: 



|P + Uki 



2eo(27r)3u;k 



F 2eo(27r)3wk 



k.(P + -ftK) 



nk 



' . P+iftK\ 1 

(E.16) 




'2eo(27r)3a;k ftA; 



sm ( 



'|P- 



(E.17) 



-^('^fc)Af |a2^.32| 



'2eo(27r)3wk 



F 2eo(2^)3a;k 



|P + iSK| 



sin 




^ ^ |P + i/iK 




1 



(E.18) 



g2g2 |p - i;iK| ± 



'2eo(27r)3wk 



1 

hk 



(-k).(P - iftK) 



(E.19) 



Combining these expressions with (E^), we obtain 

/K)a/Ip, + /(^OmIf = ^ (2^)3.;k 2/^2^^ (^|P+i;.K| " £^|P-1.K| 

(k.(P + lftK))^ . (k-(P-^^K))^ . 

|P+1;.K| --V+f.K|-3c |p_,^j^| smV-l.K| 

-2|P + i;iK|csin0|p^i^K| + 2|P - ^^K|csin0|p_i^K|) (^-20) 
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The angular integration with (k.(P+^?iK))^ will provide the ^ of the angular 
integration with |P + ^?iKp. A simplification is thus possible and we have: 

-|P + i;iK|csin0|p^i,^K| + |P - ^^K|csin0|p_i^K| 



|P + i;iK|V ^ |P-lSK|V^ 



-^IP+inKI ^|P-iriK 



2 

22 1/ 24 24 

^gC 1 / m c m c 



(E.21) 



F Expliciting the second order contribution (x ^)^^^ 
to the dressing operator (x^) 



In order to evaluate the right hand side of ( |8.21| ), we have to compute the 
(x^^)^^^ operator. Anew, we do not consider the presence of an incident 
field and the expression of (x^^)*^^"* can be deduced from ( ^.lOl) . We call 
/^^^^(— K,— X) the second order contribution. We write explicitly only 

the contribution arising from /(loooo- '^^^ creation operator C^^-* requires 
one matrix element of the coupling between the Dirac and electromagnetic 
fields. We take into account the vanishing of a matrix element when more 
than two matter indices are changed and the vanishing of two elements of 
7q,/3 to obtain: 

/y(-K,-X) (F.l) 
X < 110010^i/x'i|C(i)|110000/i/i' > /ffo'ooo 

\I \j I 111" 

X < 110001^i/x;|C(i)|110000/i/i' > /ioo 
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X < 010101^ii/xi|C*^)|110000/x;u' > /ffo'ooo 

X < 101010;Ui/x;|C(i)|110000/X;u' > /ffo'ooo (F.2) 

The values of the first order 7, determined by the Weyl's rule of correspon- 
dence, have been given previously ( B.7 - [B.l^ ). 



In order to use the form of XaB^^^'^ ( |8.28| ), we have to write (F\2) in 



the appropriate variables. We write the four terms in ( [F.2| ) as the sum of 
four contributions affected by the indices corresponding to the nature of the 
Dirac and the e. m. fields. We replace in each term the 7 by its expression 
and we take the Fourier transform with respect to the variable X. We change 
K — > — K. From our definition of Fourier transform, we have to multiply 

both sides by j2:^s " and integrate over X. That procedure provides a 

function of the variable y- order to obtain a normalisable function of P, 
that function has to be multiplied by the factor A- 



fW{2) , . _ 1 -ke / ^ f 3 

/a/3110010 i^'^J - (2vr;,)3 n \2eo{2Tif ' ^ 



nil' ... \i 



MM' MiM'i -^1 



< 110010;Ui/i;|C7«|110000^M' > ATo'oqo (F-3) 



/a/3110001 (K,P) - ^^^y^^^y d 

nil' II, 11' \' 



fifj-' mim'i A'j 

x6{p -p' -HK- hk[)e''^'^-^^—e]^' ^ .X 



"■1 



e'^C'"''(P''P) < 110001/.i/.;|C(i)|110000/./.' > AToooo(F-4) 
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mm' Mim'i 



xd{-p' -p'-hK- hk'^)e-'^'^-^^—e]^>y .X(l - e~'^'^-^) 

(F.5) 



/a/3101010 (K, P) - ^2^-^ y2i^^y Xe 
,,„' ,,' \, •' 



X(5(p + p - aK + ;iki)e-^5K.X^Le]^ .X(e*i-^ - 1) 

xe'^C'"''(-P'P) < 101010/iiM'i|C(i)|110000/xM' > /rioo(F-6) 

Trough usual manipulations, we have successsively for the action of any 
^(p, p') function: 

j (fp j dV j d^X 5(p -p' -h'K + fiki) 

= {2Tifn^- I cfp I (f'p'Sip - P - ifiK + /iki)(5(p' - P + -KK) 
xekiAi-[-Vp-Vp.]5(p,p') (F.7) 

d^p y (i^^ e^^(5(p - p' - ;iK - nk'i) 

xe-^^K.Xgj^,^^.Xe^5(p,p') 

= (27r)3fi3- / d^p f (fp'5{p - P - JfiK)(5(p' - P + \h'K + ^iki) 

% J J ^ ^ 

xek;A;-[-Vp - Vp/]5(p,p') (F.8) 
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-»PX 



p'.X 



5(P',P') 



6{p' + P + ^hK)6{p' - P + ^/iK + hk[) 



-5{p' + P + ^?iK + hk[)6ip' - P + ^ftK) 



ek;A'i-[-Vp' + Vp/]5(p',p') 
(F.9) 



d> y d^p J d^X e^^5{p + p - ftK + /iki) 



= (27r)3?i3:J y" d^p j d^p 
-6{p - P - ^ftK)(5(p + P - ^hK + hki] 



6{p-P- ^hK + Ski)5(p + P - ^SK) 



ekiAi-[-Vp + Vp]5r(p,p) 
(F.IO) 



G Computation of [x^^\x'^)^^\ 0^°^] 



In our investigation, we are interested whether that the second term in ( ^.21| ) 
may compensate the divergent contribution of the first one. We focus on 
one kind of contribution, namely the connection between /iioooo(P' P') 
itself. From (8.26), this corresponds also to the connection of J^^/^Im (k, P) 
to itself. When the argument is changed by a photon wave number, no 
divergence occurs due to the assumed behaviour of /noooo(P'P') large 
values of its arguments. We have to restrict to cases with the same value 
for the argument. The possible existence of such diagonal contributions 
can be easily realized. Indeed, the change of argument due to the photon 
in the Dirac delta functions in ( [F.7| )-( p\10| ) can be compensated by the 
change of argument arising from the creation operator in ( F.3 )-( p\^ ). The 
operator G^*^) is purely diagonal. In order that the diagonal contribution of 
the commutator provides a non vanishing result, it will be necessary that 
the gradients present in ( [F.7|) - (|F.10|) acts on the diagonal operator 

To obtain the matrix elements of X^^\x~^)^'^\ we acts with x^^'^ on both 
sides of ( [F.3| )-( pl6|) , after use of the expressions (F^ )- fUOl ). We introduce 
lower indices in expressions such as [x'''''*/^*'^^]iiooio(-^' *o keep trace of 
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the origin of each contribution: the indices 110010 are the indices of the 
intermediate state present in the creation operator. The intermediate steps 
are given only for that contribution since the other ones are computed in a 
completely similar way. 



[x(°V.T]ffio(K,P) = ^Ee/(P + 5k,P - ^K) 



h V 2eo(27r) 



x5(p-P-^K + fiki)<^(p'-P + ^K)) 



1 h 



?ek,,,.[-Vp-Vp,]<5if'^^(p',p) 



X < llOOlO/iiAi'i I C(i) 1 110000/22 Ms > /nolo (G-l) 

We replace also the function /iioooo(Pi P') its expression in terms of 
^1/.!/.' ^j^^ p^) (|g35[ ) in order that both sides of the expressions are in similar 
variables. 

Ix^'^fT^lZU^,^) = ^^Ee/(P + ^K,P - ^K) 

af) 



2eo(2vr) 



x5(p-P-^K + ftki)5(p'-P + ^K)) 
x^7^k,A,.[-Vp-Vp,]5if^^(p',p) 



X < 110010/ii/x'i|C« 1 110000/22/^2 > /'^''''H^(P-P'),^(P + P')) 

(G.2) 

The relevant first order creation operator has been determined in ap- 
pendix D and its expression ( p.l7| ) can be reported in previous formula 
( p. 21) that becomes: 

[x(°V.T]ffio(K,P) = ~Y.tt\^ + ^K,P - ^K) 

«/3 
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2eo(2vr) 



E E E / ^^^1 / I dy 



x5(p - P - + ftki)5(p' - P + ^K)) 



1 h 



X -ek,A, -[-Vp - Vp,]^^^'^^ (p', p) 



^ -Sip + ;»ki; p)/i^2i'^Hi(p - p'), ^(p + pO) 



(G.3) 

Using the relations ( |8.10| ) between the a's and the 6^s, we obtain, for ah 
contributions: 



[x(°V„T]So(K,P) = jEe/(P + ^K,P - ^K) 

a/3 

' ^ E E E/ rf'fci/rf'p / d'p' 



2eo(2vr) 



x5(p - P - + ftki)5(p' - P + ^K)) 
X -^yeic.A. •[- Vp - Vp,]^f (p', p) 



^ -5(p + hk,-p)f^^^^^^HliP - P'), liP + pO) 



(G.4) 



[x(°V„T]Si(K,P) = ^Ee/'^'CP + 5k,P - ^K) 

a/3 

' E E E/ ^'^'i / / '^V 



2eo(27r) 



x5(p - P - Jk)5(p' - P + + nk'i)) 
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1 



1 



,/i2 



•■ p -^5(p' + ^k;;p')/i'^^i'^H^(p-p'),^(p + p')) 

(G.5) 



a/9 



-I I \ I ■J 



Sip' + P + ^/iK)5(p' - P + ifiK + hk[] 

-5{p' + P + ^;iK + h]<.[)S{p' - P + ^fiK) 



1 h 



7e-klAl-[-Vp'+Vp,]5i7'^np',-p') 



;- 2^^^5(-p' - hk[;p)f'^-'>^'^{l{p - pO, ^(p + p')) 

(G.6) 



J 101010 V 



;i V 2£o(27r)3 



E E E / ^'^1 / d'P J d'P 



5{p-F- hiK + fiki)(5(p + P - hiK) 
-5{p - P - ^;iK)5(p + P - ^TlK + nki) 
X ^ ^ek,A, •[- Vp + Vf\€f'" (-P> P) 
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3 



(G.7) 



These expressions ( p.4 - |G.7 ) provide the four contributions to the matrix 



elements < l^l^'|x'-*^Hx ^)^^^ 1 1/^2 1/^2 >• particular, we have the link 
between [f^J''^^]^^^'^ (K,P) and its Fourier transform /iiqqqq- We can eval- 
uate therefore the diagonal action of [x*-'^Hx~^)*'^\ on /fi^ooo- Iiideed, 
the operator Q^^^ is nothing but Cd that is diagonal in the basis of states 
|110000;U2//2 >• In the evaluation of the commutator 0^°^], at 

least one of the operators present in X^^Hx~^)^'^^ have to act on G^''^ to pro- 
vide a non vanishing result. These operators are the substitution operators 
and the gradients Vp, Vp', Vp and Vp,. 

G.l The first contribution with the intermediate state 110010 



Let us examine more closely the contribution arising from (G.4) and compute 
the commutator with Cd- We have: 

< 1 10000^2^2 1 -^D 1 1 10000^2 /W2 >= ^ {Ep - Epi) (G.8) 



x^Ee/(P + ^K,P-^K)( 

0/3 \ 



10010= 



ih 



EEE 



h Y 2eo(27r)3 

X j d^ki J d^p J d^p'5ip - P - + hki)5{p' - P + ^K)) 
x^7eka.-[-Vp-Vp,]5if'^^(p',p) 



2eo(27r)3wk^ 
1 



i=l,3 

S{p + hki] p 



1^2 

(G.9) 
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< l^,l^,'\x^^\x-^i^^CD\llmm >nooio= -JZ^I'p'' + 1^^^^ - i^^) 

a(3 



I —I 

X 



h y 2eo(27r)2 



E E E / ^'^1 / d'p J dh' 



x6{p - P - + hk,)6{p' - P + ^K)) 
x-^^yek.A.-I-Vp - Vp,]5if ^^(p',p) 



^ -S{p + hki;p)^{Ep-Ep>) (G.IO) 



< l^lf^'lx^'Hx-' f^CDllWOOO >nooio= JE'^«^'^'(P + |k,P - |k) 

a/3 



X 



x5(p-P-^K + ?.ki)<5(p'-P + ^K)) 



xih 



^ 5(p + nki;p) (G.ll) 



-^p+riki - Ep- f^^ki 

G.2 Diagonality of the contribution for 110010 in the vari- 
ables K-P 

In order to have a possible compensation, the operator [x^^\x^^)^'^\^^^'^] 
should be diagonal in the variables K-P, as the diverging contributions of 
G(2). 

We consider first the contribution for which no gradient bears on the 
energy difference. We have to substract the energy difference (|G.8| ) for the 
initial and final states involved. The left state is for values K-P while the 
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right state involves [Ep — Epi) . After use of the substitution operator 5'(p + 
?iki;p) and the Dirac delta functions, we note that the energy difference 
vanishes identically. 

For the terms where the gradient bears on the energy difference, we need: 

V F - P"^' \7 F _ (P + ftkQc^ 

^p+ftki 

The commutator involves the computation of 

ekiAi-[-Vp- Vp,][£;p+,ki -^P'] (G.13) 



Combining ( G.ll ) and this last expression, we have the structure: 
J (fp J d^p'6{p - P - + hki)6{p' - P + ^K)) . . . 

X /iToooo(P + ^ki, pO ^ /rfo'ooo(P + ^K, P - ^K) (G.14) 

In the variables K, P, the argument of the last function is indeed K, P 
and the term could compensate the problems due to Q^"^^ that is diagonal 
in both K and P. i 

Therefore, our generalisation of Weyl's form for the quantum observables 
provide a possible compensation mechanism. We have to look at the other 
contributions arising from ( p.5 - G.7| ). 

G.3 Diagonality of the contribution for 110001 in the vari- 
ables K-P 

The contribution with the intermediate state 110001 provides a priori a 
similar contribution, with the same property regarding its diagonality in 
K-P. Indeed, instead of ( p.l4|) , we have: 

I d'p I d^p'5{p - P - ^K)5{p' - P + + Hk[)) . . . 
x/rroooo(P, P' + f^K) - /Tfoooo(P + ^K, P - ^K) (G.15) 
that is diagonal, as expected. 



permutation of the position of Axi(r) and Ax,.(r) in (5.8) would lead to a vanishing 
diagonal contribution 
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G.4 The contribution with the intermediate state 010101 
Instead of (lOTll ), we have now 

< l^^lAx^'\x-')^'^CD\llmm >oioioi= \Y.^I^^'\^ 



af3 



-K,P--K) 

2 ' 2 ^ 



2eo(2vr)2 



E E E / d% / dV / 

f I , J J J 



5(p' + p + ^riK)5(p' - p + hiK + m^:^] 

-5(p' + P + hiK + nk;)5(p' - P + ^ftK) 

1 



x5. 



a/3 



{p',-p'){-iW 



=1,3 



1 



h 1 

5(-p' - hk[;p) 



(G.16) 
(G.17) 



-E'-p'-rik; + Epf + /iWfcj 
The commutator involves the computation of 

where the gradient bears on the energy difference. 

For that term and all the terms where all the gradients act on the energy 
difference, we have the structure: 



5{p' + P + ^hK)6{p' - P + ^hK + hk[) 



-<5(p' + P + ^hK + hk[)5{p' - P + ^hK) 
x/rfoooo(-P - ^k'l, p') ^ /rroooo(P + - ;.kl, P - - ftk'i) 



-/rfoooo(P + ^K,P-::K) 



(G.18) 



2 ' 2 

In the variables K, P, the argument of the last function is K, P and the 
term can compensate the problems due to G*^^-* that is diagonal in both K 
and P. 
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G.5 The contribution with the intermediate state 101010 

It is easily realised that we obtain a contribution diagonal in the variables 
K, P and the contribution when no gradient acts on the energy difference 
vanishes identically. Indeed, instead of ( G.18 ), we obtain 



-5{p 



6{p-P - -hK + hki)6{p 
P - ^hK)6{p + P - ^hK + hki] 



-hK) 
2 ^ 



x...5(-p-ftki;p')/So(P.P' 



d^p' J d^p'5{p 



110000 V 

P - hiK + Ski)(5(p + P - ]^KK) . . . 



110000 IP' 



-p-ftki)-/r 



ioooo(P+2K-^ki,P--K 



-/rroooo(P + 5K,P-::K) 



2 



(G.19) 



H Contribution to the 0^ 

We now turn to the explicit evaluation of the diagonal singular part of 
(^(0)$(2)(^-i){0) ^0(2) ^ [x(0)(x-i)(2),e(0)]. The diagonality has to be 
considered with respect to three criterions. In the first one, we consider 
in 0^^) and the diagonal connection between the state 110000 and 

itself. We expect indeed that eventual divergences in that connection are 
independent of the other contributions. In the second one, we consider the 
contribution diagonal with respect to the value of the argument of the func- 
tions. The third criterion refers to the spin variables. We will not compute 
the matrix elements for all spin indices but introduce usual functions describ- 
ing the spin properties of one particle. It is useful at this stage to introduce 
the distribution function /(K,P,t), describing the electron independently 
of its spin orientation: 

/(K,P,t)=^/ffoooo(K,P,t) (H.l) 

(T 

and the spin orientation distribution function (in short SODF) (more pre- 
cisely, it is a vector but we keep the terminology of the litterature |Q , [^] 
) A4(K, P, t) defined in terms of the Pauli matrices trough a partial trace 
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over the spin variables: 



M{K,P,t) = tropins a(/iioooo(K,P,t)) (H.2) 

where /iioooo(K, P, t) is the one electron Wdf considered as a matrix in spin 
variables. 

These relations ( p?ll|L2| ) can be inverted to provide: 

ATo'ooo(K,P,t) = \[f(K,P,t)I+M(K,P,t).aU, (H.3) 
The diagonal kinetic operator acting on this new set of functions receives 

(2) 

a M index. The operator (p\/[ is the operator acting on those functions 
/(K, P, t) that contributes to their time derivative. Vectorial (in spin space) 
and tensorial operators would also be required to complete the kinetic equa- 
tions, that involve also A<(K, P, t). Their explicit form is irrelevant here. 
We focus here on the diagonal contributions arising from [x^^\x~^)^'^\ Q^'^^]. 
We have the elements to compute the diagonal contributions leading from 
110000 to itself. We place an index a to notify the fact that only that con- 
tribution is included and we still note the intermediary state explicitly. For 
the first contribution, the required expressions in the K and P variables can 
be obtained from ( |GHD , (|G1| ) and ( |Gl| ): 

[x(°)(x-^)(^),eW]M.nooio= E ^^E'^S^P + 5k,P-^K) 



\^ h Y 2eo(27r)3y ^7 ^/uJk^ i 

y-Qec-^ E (et, A, + |k - Bki.p + |k) 

j=l,3 

^ (H.4) 



The other contributions are determied in a similar way. From ( p.q ), we 
obtain: 

[x(°)(x-^)(^\eW]M.noooi= E ^^E^^";"(P + 5k,P-^K) 
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d^k[ — ii-ek/ v^. -Vp (-Bp_^|K - ^P-^k] 



^^iA.2iK(p _ ^K,P - -K - ^ik'i)— (H.5) 

2 2 --E^P-|K + -^P-fK-sk; + ^^fc'i 

From ( p. 61) (second contribution with the 66, hence a (—1) factor): 



- ^p iK 



x# ''^^(P + iftK + ;.k;, P + IhK)- -— ^ —— (H.6) 

From ( p. 71) , (second contribution with the 66, hence a (—1) factor): 

[x(°)(x-^)(^\eW]M.ioioio= E ^^E^«';^p + 5k,p-^k) 



xjjf "'^'^P - ^ftK,P - iftK + nki)— -i ^ (H.7) 

Elementary but lengthy algebra provides the relations (v is an arbitrary unit 
vector, Vs its s component) 

E E E v.E'^«'/'(P'P')'^«"/'"np',P")4f ''^^(p",p) 

= 8sin^cos^(p.v) + 8sin^cos^(p".v) (H.8) 
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E E E vsE'^a";"(p,p')e/'"^(p'>p")4"^'"^(p">p) 

a ■ dp dp - ^ o ■ dp" , . „ . . 

= 8sin ^ COS ^(p.v) — 8sm -|- cos -|-(p .v) 



(H.9) 



For the expression ( [H.^ ), we have to consider v = ej^^;^^. We have also 



P+fK 



Vp-E/p_ 



P — -K 



P-|K 



(H.IO) 



We deduce for (|tL4|): 

A/allOOlO 



sm ■ 



V2eo(2 

P+|K| V+fKI /P+|K 



cos ■ 



|P + |K| 



. f^lP+^K-fikil 
- sm = cos ■ 



'lP+§K-riki| / P + |K-nki 



|P + |K - nki 



xekiAi- 



„P + |K oP-|K' 

-C^^; ^ h C^^; ^ 



E- 



P+fK 



^P |K 



-^P+f K ~ -^P+f K-;iki ~ ^^fci 

(H.U) 



and the summation over the polarisation index Ai can be performed easily, 
using (only the transverse components play a role): 

E (a-ekiAi)(b.ek^Ai) = (a.b) - (a.ki)(ki.b) (H.12) 

Ai=e,e' 



Afal 10010 




\ 2eo(27r 



COS 




P + |K _^ P - |K' 



P+fK 



E. 



P + |K 



E- 



+ 



2 
P 



K 



P+fK 



E. 
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^|P+^K-;iki| ^|P+fK-Rki| 
+ sin cos 



P + |K 



|P + |K 



P + |K ^ P - |K 



rK 



P + |K 



P+|K P-^K 



E- 



+ 



P+|K 



E- 



P-|K 



.ki. 



^P+fK 



(H.13) 



For large uJki , taking into account that angular integrations may vanish, the 
integrand behaves as 



X sm ■ 



MallOOlOlrfom 



■ cos ■ 



^|P+fK 



2eo(27r)3 

h 



P + fK 

|P + |K| 



P + |K _^ P - |K' 



E- 



^P-|K 



(H.14) 



From ( p. 61) , we obtain: 

A/aOlOlOl 




V2eo(2vr 
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E- 



(H.15) 



and 



-^('^fcjA/aOlOlOl I rfom 



( 



1 2 



y2eo(27r)3y cj^^ 2hujk^ 3 



^1 



X sm ■ 



P+#K| 



■ COS ■ 



V+fKI / P + |K 



|P + |K| 



P+|K P-|K 
-^P+fK ^P-|K 



(H.16) 



The linear divergence cancels out from the two expressions ( |H.14[ ) and 

(MM)- 

From and (|HJ), we obtain, using the hermiticity of the 6 matrices: 

s=l,3Ati:/^2 



8sm — COS — (p.vj + osm COS -^(p .vj 
2 2 2 2 



(H.17) 



E E E v,^e;"(p'.p)^^^'"^(p">p')4^^'^^(P,P": 

a,/3 s=l,3/ii,/i2 M 



ism COS — (p.vj — 8sin COS -|-(p .vj 



(H.18) 



In (|G.6| ), we can use ( H.17 ) with the correspondence p' — > P + ^K, p 



P - |K, p" ^ P - |K - hk[. Therefore, 



[^(O)(^-l)(2)^0(O)] 



A/allOOOl 



\2eo{2 



TT 



I . '^|P-|K| 
X < sm — cos ■ 



^|P-fK| 



P - |K 
JP-|K| 

, ''iP-fKl 



P + |K _^ P - |K' 



E- 



P+fK 



E- 



P-|K 



P + |K _^ P- |K' 



^P+fK 



^1 



P-fK 



^|P-4K+;ik'| ^|P--|K+?ik' 
+ sm cos 
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P-|K 



-K 



hk[\ 



P + |K _^ P - |K' 



P+|K 



-K 



P-|K 




-K 

2^^ 



Et 



P-|K 



+ 



-K 



E- 



P-|K 



^p-f K + ^p-|K-nk; + ^^fc; 

H.13|) by replacing K 



This expression is obtained from 



The bracket 



P+fK _^ P- 



(H.19) 



-K, ki ^k'l. 



generates a change of sign, analogous to 



the change of sign required in the correspondence in the elements of 0^^). 
In ( p.?! ), we can use ( |H.l^ ) with the correspondence p' ^ P + |K, 



|K, p" 



|K 



fiki. We then obtain 



[^(O)(^-l)(2)^0(O)] 



A/alOlOlO 




X < sm ■ 



|P-iK| i^lP- 
COS 



rKI 



P-|K 
|P-|K| 



ki.: 



P + |K P-|K 



P+fK 



E- 



P-|K 




sm ■ 



P-|K+;iki| 



P + |K P 

— — + 



-K 

2^ 



E, 







■ cos ■ 



^P+fK 
IP-fK+nki 



P-fK 



P - |K + ftki 
|P-|K + ftki| 



P + |K _^ P - |K' 



E. 



ki-T 



K + hki 



fK + ;iki| 



+ tK 
P + §K 



+ 



2 
P 



K 



-K 

2^ 



E- 



P |K 



•ki 



-E- 



(H.20) 



For the dominant part, linearly divergent, we have the same cancellation for 
the sum of gT9| ) and ( gL20|) as for (^T|) and ( gLl^ ). 

Therefore, the remaining contribution diverges logarithmically. 
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